United States Patent 

Moose 



US005 166924 A 

[ii] Patent Number: 
[45] Date of Patent: 



5,166,924 
Nov. 24, 1992 



[54] ECHO CANCELLATION IN 

MULTI-FREQUENCY DIFFERENTIALLY 
ENCODED DIGITAL COMMUNICATIONS 

[75] Inventor: Paul H. Moose, Carmel, Calif. 

[73] Assignee: Mercury Digital Communications, 
Inc., Monterey, Calif. 

[21] Appl. No.: 566,290 

[22] Filed: Aug. 10, 1990 



Related U.S. Application Data 

[63] Continuation-in-part of Ser. No. 490,769, Mar. 6, 1 990, 
Pat. No. 5,063,574. 

[51] Int. CI.' H04L 5/14; H04M 9/08 

[52] U.S. CI 370/32.1; 379/410 

[58] Field of Search 370/32.1, 32; 379/410, 

379/411 



[56] References Cited 

U.S. PATENT DOCUMENTS 



4,799,214 
4,868,874 



1/1989 
9/1989 



Kaku 370/32.1 

Takatori et al 370/32.1 



Primary Examiner— Douglas W. Olms 
Assistant Examiner— Wellington Chin 
Attorney, Agent, or Firm— Davis & Schroeder 

[57] ABSTRACT 

A differentially encoded digital signal waveform is 
generated as a discrete time representation of a desired 
analog signal utilizing multi-frequency modulation 
techniques. The computational capability of present 
day, industry-standard microcomputers equipped with 
a floating point array processor or digital signal proces- 
sor chip is utilized to perform digital frequency encod- 
ing and compute both discrete Fourier transforms and 
inverse discrete Fourier transforms to provide a trans- 
mitter and receiver system utilizing suitably pro- 
grammed microcomputers coupled by a communica- 
tions channel. 
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combination of phase and amplitude modulation known 
ECHO CANCELLATION IN MULTI-FREQUENCY as quadrature amplitude modulation (QAM). For exam- 
DIFFERENTIALLY ENCODED DIGITAL pie, transmission of data over telephone lines designed 

COMMUNICATIONS to carry analog voice frequency signals (VF) is re- 

5 stricted to a frequency band of about 300 to 3500 Hz. 
CROSS REFERENCE TO RELATED Digital modulation of an 1800 Hz carrier frequency 

APPLICATION using QPSK (2 bits per symbol) or 16-QAM (4 bits per 

This application is a continuation-in-part of copend- symbol) utilizing a modulator-demodulator device 
ing application Ser. No. 07/490,769 filed on Mar. 6, (modem) provides data communication at 4800 and 
1990, now U. S. Pat. No. 5,063,574. 10 9600 bits/sec, respectively, when transmitting at a rate 

of 2400 symbols per second. The unprocessed modula- 
BACKGROUND OF THE INVENTION tor ^ct signal thus provided has a bandwidth sub- 

The present invention relates generally to communi- stantially greater than the available channel bandwidth 
cations of digital signals over bandlimited, unequalized, so it must be filtered to fit in the available band. These 
noisy transmission channels at high data rates and, more 15 analog filters must be carefully designed so as not to 
particularly to the use of microcomputers such as per- introduce smearing, known as inter-symbol interference 
sonal computers and the like to accomplish multi-fre- (ISI) between adjacent symbol waveforms or the sym- 
quency modulation and demodulation for transmission bols will be decoded in error. A raised cosine filter 
of a multitude of frequencies over a single channel. function is typically employed. Even then, the filtering 

In the modern world, computer to computer is 20 action of the bandlimited VF channel will introduce ISI 
quickly becoming the fundamental link for message and <j ue to non-linearities (group delay) in its phase re- 
data traffic; i.e., information exchange and propagation. sponse. These non-linearities are most pronounced at 
In some cases the link is one of many in an extensive t fc e band edges so that the symbol waveforms must be 
network; in others, it may be represented by a single filtered to the 2400 Hz band from 600 to 3000 Hz before 
point-to-point communications link. The physical me- 25 ^ Qver the channel The receiving section of 

dium constituting the link for transmitting the data may ^ demodulator contains an adaptive filter known as an 
be wire or optical fiber. Or, it may be a microwave or equalizer ^ orc jer to remove any residual ISI intro- 
satellitelink utilizing radio frequency propagation. (Sat- duccd by the dclay ^ the ^ t0 3000 Hz band, 
ellites are also being developed for broadcasting digital ^ equa]izer must ^ ^tined to the particular group 
audio and digital television signals^ In all cases, the 30 toychlBCtetorf ^ 
actual signal used .to carry the data bits must be a prop- ^ he{o data can * 1™^^. 
erly modulated analog signal with sufficient energy and whkh lization ^ rcduce the 

the objective is to ,rans^ information at a high rate 35 ^^^^ 
with a low rate of errors from the transmitter to the 
receiver, 



receiver and hence limits the number of bits that can be 
encoded into each baud. A baud is a digitally encoded 



Prior art modulation methods for bandpass channels «y«*d wavef °™; In « ha ? n found that the 

utilize amplitude and/or phase modulation to carry combination of ISI and additi ve noise limit transmission 
signal information on a carrier wave in a channel fre- 40 to « ther two ' or under ldeal conditions four, bits per 
quency band. When the information source is a finite baud when transmitting 2400 bauds per second over 
alphabet, numbers, data, etc., or quantized analog 2-wire switched telephone circuits, 
sources such as digitized audio, television, or facsimile, In order t0 increase the rate at which data is transmit- 
then only a finite number of signal states are required to ^e modems employ multi-frequency modulation 
represent, or code the source. 45 (MFM) techniques. MFM utilizes multiple carrier fre- 
Digital data is frequently communicated using pulse quencies within the available bandwidth, each fre- 
code modulation techniques at baseband, i.e., in the quency independently modulated with digital lnforma- 
frequency spectrum from 0 that is to some maximum, or ™ phase and/or amplitude. The frequencies are 
upper frequency limit. If a bandpass transmission is linearly combined and transmitted as a single digitally 
required, i.e., is transmission in a frequency band be- 50 encoded waveform, termed a baud, during a finite time 
tween an upper frequency and a lower frequency, single interval called the baud interval. U.S. Pat No. 4,731,816 
side-band (SSB) modulation into and demodulation issued to Dirk Hughes-Hartogs on Mar. 15, 1988 and 
from the desired frequency band may be employed. For U.S. Pat No. 4,601,045 issued to Daniel P. Lubarsky on 
example, polar voltage pulses at 48K bits/sec are fil- Jul. 15, 1986 disclose examples of modems employing 
tered to attenuate frequencies above 36 KHz and SSB 55 MFM techniques. Hughes-Hartogs teaches mininiizing 
modulated onto a 100 KHz carrier to fit into the band- inter-baud interference by introducing a small guard 
pass channel from 60-104 KHz. time between successive bauds during which no signal 
Frequency modulation techniques are also employed is sent to prevent received baud overlap. Lubarsky 
for digital data. Typically, these involve sending one of teaches concentrating most of the signal in the center of 
two, frequency shift keying (FSK), or one of M (M-ary 60 the baud interval thus causing the signal to taper off to 
FSK), frequencies spaced across the available fre- zero near the ends of the baud and minimizing interfer- 
quency band and may be used in applications where ence between baud intervals by eliminating abrupt 
bandwidth efficiency is unimportant as they operate at changes. If the baud interval is of sufficient duration 
less than one bit per Hz of available bandwidth. compared to the guard time or taper time, then the loss 

Alternatively, digital information may be encoded 65 of data rate is relatively insignificant, 
onto a basically analog carrier frequency, centered in If the bauds are to be long in duration, then many bits, 
the available frequency band, using phase modulation and consequently many frequencies, may be transmitted 
(PSK), differential phase modulation (DPSK), or a in each baud. However, in order to prevent inter-fre- 
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quency interference, the frequencies must not be spaced tions channels. The present invention functions in pack- 
too closely. Ideally, the frequencies are made orthogo- et-switched, circuit-switched, or broadcast enviorn- 
nal over one baud interval. This is known as orthogonal ment and operates with either baseband or bandpass 
frequency division multiplexing (OFDM). A frequency channels. An MFM digital signal waveform is gener- 
set will be mutually orthogonal if the frequencies are 5 ated as a discrete time, i.e., sampled data, representation 
separated at multiples of the reciprocal of the baud of a desired MFM signal using an industry standard 
interval. For example, a system employing a rate of 10 microcomputer, such as a personal computer or the like, 
bauds per second in the VF band from 300 to 3500 Hz equipped with a floating point array processor or a 
could contain 320 frequencies spaced 10 Hz apart. If digital signal processing (DSP) chip and utilizing a 
each frequency is encoded with 4 bits of information, 10 computer software program. For signal transmission, 
then 1280 bits will be encoded into each baud and the the discrete time signal samples, which are stored in 
system will have a throughput rate of 12,800 bits/sec. random access computer memory (RAM) as digital 

The major disadvantage associated with the MFM words in a time domain array (TDA), are clocked from 
(OFDM) systems described above is that although they memory via a suitable input/output (I/O) port of the 
can effectively eliminate inter-baud interference and 15 microcomputer at a clock rate which is a large multiple 
inter-frequency interference within a baud, they must of the baud rate and which is at least twice the highest 
be demodulated using fully coherent receivers for each frequency to be sent in the MFM analog signal. The 
frequency in order to obtain the high data rates desired. digital words are coupled to a digital-to-analog con- 
Since the multitude of frequencies are subject to differ- verter (DAC) as a serial data stream for conversion to 
ent and unknown amplitude and phase changes intro- 20 the desired analog MFM signal for transmission, 
duced by the transmission channel, these channel prop- The computer program creates the discrete time sig- 
erties must be measured during the initiation phase of nal as follows. Digital data to be transmitted are ac- 
the communication process and prior to transmission of cepted from a message file n bits at a time, differentially 
data. One example of this technique is disclosed by encoded and converted to real and imaginary signal 
Hughes-Hartogs. Such methods are exceedingly com- 25 values in accordance with the chosen signal modulation 
plex and computationally intensive easily rivaling the constallation for each of the tones. Each encoded n bit 
cost and complexity of adaptive equalization. Further- work is stored in a position in the lower half of a corn- 
more, such OFDM systems will not be effective for plex valued frequency domain array (FDA). The corn- 
direct MFM signalling in bandpass systems such as the plex conjugate of the encoded value is stored in its 
60-104 KHz example described above or in a prototype 30 image position in the upper half of the FDA. Each 
model of a UHF satellite sound broadcasting system as position in the array is assigned a position number 
discussed by Alard et al in "A New System of Sound which corresponds to a harmonic of the baud rate of the 
Broadcasting to Mobile Receivers", presented at the system. Only those array positions corresponding to 
"Centre Commun d 'Etude de Telecommunication et desired tones in the MFM signal are filled with modula- 
Telediffusion" in France and published by the IEEE in 35 tion values. All other array position values are set to 
1988. zero. 

Differential encoding of the carrier frequencies pro- In the preferred embodiment, K encoded words are 

vides a practical solution to this problem with an atten- placed in adjacent array positions beginning at position 

dant 3-db loss of signal-to-noise ratio performance kj + 1 and ending at position ki 4* K. Reference values 
against additive noise. The conventional method to 40 for differential decoding are placed in position ki; this 

differentially encode information is from baud to baud results in an analog signal strictly bandlimited to the 

as is customary in conventional DPSK and as done in band between ki — 1 times the baud rate and ki+K+1 

the OFDM system disclosed by Alard et al. However, times the baud rate. For example, a baud rate of 10 per 

this method produces undesirable results when the baud second, a kj of 31 and a K of 320 results in a modulated 

interval is long, due to channel instability such as that 45 analog signal strictly bandlimited to the band between 

introduced by fading. Further if asynchronous or 300 Hz and 3520 Hz. A baud rate of 250 per second, ki 

packet transmissions are utilized there may be a signifi- of 241 and K of 174 results in a modulated analog signal 

cant reduction of data rate when only two or three strictly bandlimited to the band between 60,000 Hz and 

bauds are sent since differential encoding in time re- 104,000 Hz. A baud rate of 30 per second, a kj of 537 
quires utilizing one baud as a reference. 50 and a K of 128 results in a modulated analog signal 

Frequency differential encoding of multiple carriers strictly bandlimited to the band between 16,110 Hz and 

was utilized in the HF modem disclosed by Gene C. 19,980 Hz. 

Porter in "Error Distribution and Diversity Perfor- A real valued TDA is obtained by operating on the 

mance of a Frequency-Differential PSK HF Modem", FDA with an inverse discrete Fourier transform 

IEEE TRANS. ON COMMUNICATION TECH- 55 (IDFT) algorithm. The values in the TDA are guaran- 

NOLOGY, 16-4, August 1968, pages 567-575. While teed to be real because of the complex conjugate image 

frequency differential encoding minimized many of the symmetry established in the FDA. The total number of 

above described problems encountered with OFDM real values in the TDA, designated as k x , is the same as 

systems, the circuitry required for generation and de- the total number of complex values in the FDA. In a 
modulation of the signals was unduly complicated and 60 preferred embodiment, k* is made a power of two and 

did not reliably maintain the necessary orthogonality the IDFT is executed using a fast Fourier transform 

between the carrier frequencies. (FFT) algorithm. The TDA contains the discrete time 

CTT ww APV nP TTLip rwPMTinw "S™ 1 representation of one baud of digitally encoded 

SUMMARY OF THE INVENTION flata ^ ^ of ^ superposition of kx/2- 1 tones, 

The present invention provides a method and appara- 65 or carrier frequencies, each modulated with a phase and 

tus utilizing multi-frequency modulation (MFM) tech- amplitude corresponding to the polar representation of 

niques for modulating and demodulating digital infor- the complex number that was stored at its harmonic 

mation signals for use with bandlimited communica- number in the FDA as described above. The tones are 
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harmonics of the baud rate and are mutually orthogonal with more points are selected, a guard time equal to the 

on the baud interval. The TDA values are clocked out group delay with zero signal values is inserted between 

through the microcomputer I/O port to the DAC at a bauds in the TDA. 

clock rate of k x times the baud rate. Demodulation of MFM communication signals is 
In the preferred embodiment where an adjacent sub- 5 accomplished in a manner that is the inverse of their 
set of K of the k x /2 — 1 available tones are modulated generation. The received analog signal is sampled by a 
with information and the others given amplitudes of clock having the same frequency as the transmit clock, 
zero, nK bits are transmitted in each baud. In the 10 Sample values are analog-to-digital (ADC) converted 
baud per second example containing 320 adjacent tones and stored in the receive personal computer memory in 
between 300 Hz and 3520 Hz, cited above, k x of 1024 10 a receive TDA. The data are extracted from the TDP 
(2 10 ) leads to a clock rate of 10,240 Hz. With n=4 bits one baud at a time. The receive computer program 
encoded into each tone, 1280 bits are sent with each executes a k x point discrete Fourier transform (DFT) 
baud for a data rate of 12,800 bits/sec. In the 250 baud creating a receive FDA. Complex values from the posi- 
per second example containing 174 adjacent tones be- tion values in the FDA corresponding to the K+l 
tween 60 KHz and 104 KHz tones, k, of 1024 leads to a 15 tones in the MFM signal are differentially decoded 
clock rate of 256 KHz. With n=2 bits encoded with producing K complex values that are stored in a receive 
each tone, 348 bits of information are sent with each digital signal array (DSA). These values are demodu- 
baud for a data rate of 87,000 bits/sec. In the 30 baud lated in accordance with the detection algorithm of the 
per second example containing 128 tones in the band of selected modulation constellation to produce the K, n 
1 6, 1 10 Hz to 1 9,980 Hz a k x of 2048 leads to a clock rate 20 bit words of the baud. The entire packet is decoded by 
of 61,440 Hz. With n=3 bits per tone, 384 bits of infor- repeating the process for each of the L bauds, 
mation are transmitted with each baud for a data rate of One advantage of differential encoding in frequency 
11,520 bits/sec. is that the receiver clock frequency need not be phase 
It should be noted from the examples above, and it is coherent with the transmit clock frequency. Further, 
a significant feature of the present invention, that band- 25 small frequency offset errors, on the order of one part in 
pass and baseband signals are generated in an identical 10 s , are acceptable with no degradation in performance, 
manner; no further modulation steps are required. The first baud in each packet is utilized as a synchroni- 
ln the preferred embodiment, a time domain packet zation baud to establish baud synchronization. Because 
or array (TDA) is generated for transmission by the MFM signals have the statistical characteristics of band- 
computer program of the transmitting microcomputer 30 pass white noise, the autocorrelation function of an 
consisting of the TDA of a synchronization baud fol- MFM signal with random input data has a very sharp 
lowed by the TDAs of L data bauds. A packet transmits spike for its central peak. Thus, a polarity-only matched 
L bauds of K tones each modulated with n bits per baud filter matched to a pre-selected synchronization baud 
or nLK bits per packet at a data rate of n bits/Hz of which generates a sharp spike at the end of the synchro- 
occupied channel bandwidth. Overhead comprises one 35 nization baud may be used to initiate data transfer from 
baud for synchronization and one tone for differentially the ADC to the receive TDP via the selected receive 
encoding the data in frequency. microcomputer computer I/O port. In a preferred em- 
Frequency domain differential encoding eliminates bodiment, data is transferred via the direct memory 
the need for correcting the channel amplitude and phase access (DMA) channels in both the transmit and receive 
response, i.e., frequency domain equalization as re- 40 microcomputers using the byte transfer mode. For very 
quired by the prior art, in the following manner. The high speed asynchronous operation, block DMA trans- 
tones are spaced at the baud rate. By using a very low fers are preferred. 

baud rate, the tones are very close together and there- In some data communications settings, such as corn- 
fore experience substantially the same group delay and puter-to-computer links, full duplex transmission is re- 
attenuation when transmitted through the channel. In 45 quired, i.e., information nodes which transmit and re- 
the present invention, phase (and/or amplitude) is dif- ceive simultaneously. But when operating in the full- 
ferentially encoded from tone-to-tone. Differential de- duplex mode over the same physical communications 
modulation at the receiver cancels out unknown phase medium using the same channel bandwidth and time 
(and/or attenuation) introduced by the channel. Bulk slots, a node's own transmission can interfere with its 
delay through the channel is accounted for at the re- 50 reception from the other node. This is generally known 
ceiver with the synchronization baud using a polarity as listener echo and is a particular problem in dial-up 
only matched filter to establish the timing reference for two-wire public switched telephone circuits. In order to 
the received baud signal samples. combat this problem, high-speed dial-up modems are 
One advantage of utilizing a very low baud rate, i.e., equipped with echo cancellers, devices that remove that 
long baud intervals, is that the need for a guard time 55 component of a node's received signal that is due to its 
between bauds, or for tapering the baud signal ampli- own transmission. As much as 75 db of echo suppres- 
tudes at their beginning and end, i.e., "shoring" as dis- sion may be required in order for a system to achieve 8 
closed by Lubarsky; is reduced or eliminated. This is bits per Hz of bandwidth efficiency at a signal-to-echo 
because only a very small fraction of the baud energy is ratio of —27 db. The present invention describes a dual- 
involved in the inter-baud interference. For example, 60 stage, data driven echo cancellation algorithm for 
over unequalized telephone lines the variation in delay MFM that achieves up to 75 db of rejection of near-end 
for tones from 300 Hz to 3500 is only several millisec- echo and up to 35 db of rejection of far-end echo, 
onds. For bauds of 100 milliseconds or more in length The present invention provides an MFM transmitter 
only a small percentage of the baud energy is involved and receiver which utilizes the computing power avail- 
in inter-baud interference. This is not significant when 65 able in state-of-the-art industrial standard microcomput- 
demodulating QPSK constellations (2 bits encoded per ers to encode, modulate, demodulate and decode digital 
tone). Under particularly severe unequalized group data utilizing discrete Fourier transform (DFT) tech- 
delay conditions and or when modulation constellations niques. In a data communications system comprising a 
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transmitting microcomputer and a receive microcom- FIG. 13 is a flow diagram illustrating demodulation 
puter linked together via a communications channel, an and decoding for multi-frequency differential quadra- 
inverse discrete Fourier transform (IDFT) m the trans- ture phase shift keyed modulation in accordance with 
mit computer and a DFT in the rec eiver c omputer are tne principles of the present invention; 
executed as Fast Fourier Transforms (FFT) in real time, 5 FIG> 14a ^ a illustrating the transition states 
the only requirement being that each microcomputer be for mu i t i fc f re quency differential l^quadrature ampli- 
equipped with readily available floating point array tU(Je modulatjon to th e principles of the pres- 
processors. Each microcomputer is coupled to the com- ^ mvent j on . 

munications channel via signal converters such as „ T — -.- t .* .„ * * * t.™ 

readily available ADC and DAC, respectively, chips. 10 J 10 '™" " '*f gn ? f^^g ^ decoded phase 

Thus, the present invention provides a "universal J«« ;** differential encodmg strategy shown m 

modem" in which the bandwidth and data rate are eas- FIGl . . . 

ily controlled by changing parameters in transmit and FIG - 15 ls a diagram illustrating a decoding decision 

receive software programs. trce f° r decoding differentially encoded data in accor- 

Differential encoding of the phase from tone-to-tone 15 dance with the encoding strategy shown in FIGS. 14a 

eliminates the requirement for a phase-locked or coher- and 14b; 

ent clock reference between the transmit and receive FIG. 16 is a flow chart illustrating a transmit algo- 

terminals and greatly reduces the signal sensitivity to rithm for use in the transmit microcomputer shown in 

phase variations introduced in the physical portions of FIG. 5; 

the communications link. Baud synchronization is ob- 20 FIG. 17 is a flow chart illustrating another transmit 

tained at the receiving terminal utilizing a simple polari- algorithm for use in the transmit microcomputer shown 

ty-only correlator circuit operating on a synchroniza- m pjQ^ 5. ^ 

tion baud transmitted as the first baud in each data FIG. 18 is a flow chart illustrating a receive algo- 
packet. The amplitudes of the individual tones in the rithm ^ m the receive ^0^^ shown in FIG. 
MFM signal can be controlled at the transmit FDA to 25 5 
provide equalization of the frequency-dependent signal- 
to-noise ratio which may be present in frequency selec- DETAILED DESCRIPTION OF THE 
tive channels. PREFERRED EMBODIMENT 

BRIEF DESCRIPTION OF THE DRAWING 3Q MULTI-FREQUENCY MODULATION 

A fuller understanding of the present invention will Referring now to FIG. 1, an analog multi-frequency ] 

become apparent from the following detailed descrip- modulation (MFM) transmit signal generated according 

tion taken in conjunction with the accompanying draw- t o the principles of the present invention comprises 

ing which forms a part of the specification and in which: "packets" of multiple tones or frequencies which are 

FIG. 1 is a diagram illustrating the relationship be- 35 differentially modulated in amplitude and/or phase : 

tW ^^? e . qUenC Af d time f ° r ™ signal packet; adjacent frequenc jes. The tones 13 are present 

FIGS. 2a and 2b are conceptual block diagrams illus- simu)taneously durin ^ a submter val of the packet 11 

rating the transmitter and receiver sections, respec M fl baud 15 The kets u ^ ^ localed : 

tively, for an MFM communication system according . . „ . . . v 4 . . . 

to the principles of the present invention; 40 "^"J* ™ the Kf^"* 10 tune ' 

FIG. 3 is a diagram illustrating the complex envelope W,th furthcr reference to FIG. 1, the following terms 

for one QPSK tone; are defined: 

FIGS. 4a-4d are phase plot diagrams illustrating the T: Packet length in seconds 

phase encoding for DQPSK encoding; AT: Baud length in seconds 

FIG. 5 is a functional block diagram of an MFM 45 Baud length in number of samples 

communications system according to the principles of L; Number of bauds per packet 

the present invention; At: Time between samples in seconds 

FIG. 6 is a functional block diagram of the transmit- fjt=l/At: Sampling or clock frequency for digjtal-to- 

ter block shown in the block diagram of FIGS. 2a and analog and analog-to-digital conversion in Hz. 

5; 50 Af=l/AT: Frequency spacing (minimum) between 

FIG. 7 is a functional block diagram of the receiver MFM tones, 

block shown in the system of FIGS. 2b and 5; K: Number of MFM tones. 

FIG. 8 is a functional block diagram of the data ac- Since At = ATA* the sample frequency f*=k x Af. 

quisition board shown in FIG. 7; Consequently, there are a maximum of kx/2-1 tones 

FIG. 9 is a functional block diagram of a data-driven 55 aced at mtcrva]s of Af m Af a md lcss 

echo cancellation circuit for a full-duplex communica- than { n m ^ N ^ f ■ ^ ^ 

pon system according to the principles of the present de ^ ^ ormati ^ du ^ g ^ Uu 7£ 

invention; „ - - A . . , 

FIG. 10 is a diagram illustrating the first stage echo ^^^M^^^^^^^^ 

cancellation training algorithm for the echo cancella- 60 ta ™ an iP lltudes °f *»> dunn « any or all bauds 

tion circuit shown in FIG. 9; of ±c P"**- For example, to generate bandpass signals 

FIG. 11 is a diagram illustrating generation of the between frequencies fi and ft, only tones between har- 
weighting coefficients for the second-stage ttaining monics ki«fi/Af and k 2 =ft/Af will be allowed non- 
algorithm for the echo cancellation circuit shown in zero amplitudes. Here, the m axi m u m number of tones 
FIG. 9; 65 will be K=k2— kj + 1 and the signal bandwidth will be 

FIG. 12 is a flow diagram illustrating multi-frequency W= KAf. 

differential quadrature phase shift keyed modulation in Mathematically the \ th baud of the analog transmit 

accordance with the principles of the present invention; signal can be described by; 
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x/u) - 2 xidu) 



0) 



where, 

x\i<{u)=A\kCO& (2irkAfu+4> lk )$£u£T. (2) 

In equations (I) and (2), u is time referenced to the 
beginning of the baud 15. Actual real time is 
t=to+ 1AT+ u where to is the time of initiation of the 
0 th baud, i.e., the beginning time of the "packet". The 
discrete time signal corresponding to the \ th baud 15 is 
found by sampling equations (1) and (2) at the sampling 
intervals At= l/f x and is given by; 



10 



15 



TABLE I-continued 


k 


Re{X(k)} 


Im{X(k)} 


l 


0 


0 


2 


0 


0 


3 


XR3 


XI 3 


4 


XR4 


XI4 


5 ' 


XR5 


X15 


6 


0 


0 


7 


0 


0 


8 


0 


0 


9 


0 


0 


10 


0 


0 


11 


XR5 


-XI5 


12 


XR4 


-XI4 


13 


XR3 


-X13 


14 


0 


0 


15 


0 


0 



x&) = 2 xifc) (3) Referring now also to FIGS. 2a and 2b t a block dia- 

* gram of a transmitter section 21 and a receiver section 

20 23, respectively, are shown. The k x point DFT of xi(n), 
where, as defmed by equations (5) and (6), i.e., the complex 

valued FDA is loaded on line 22 into the IDFT block 25 
*i*<«)=^ cos (2ir*«/*,+*uW)SiiS*,-i. (4) tQ provide a ^ of discre te time samples in the TDA 

. J . r _ corresponding to the desired MFM analog transmit 

In equations (3) and (4), n is discrete time referenced 25 d ^ ^ ]es ^ clockcd out ^ DA C 27 at a 

to the beginning of the baud 15 In general, k may take ^ f u of & %Q ide ^ desired ^ 

on all integer values between 1 W2- 1; k is re- transmit si j Xl(u) on line M A lete ^ 

ferredtoasthe'W^ £ ^ cd b L titions J thc a £ ve 

frequency kAf A baud 15 interval he., a tune AT, con- ^ * $ 

tains exactly k cycles of tone k. Thus adjacent : tones 30 Demodulation of MFM is the inverse of the signal 

differ by one m the number of cycles they make during ^ recejved ^ ^ 

8 The k/point Discrete Fourier Transform (DFT) of ? n lme 36 * the re f eiver 23 is ^P 1 * at * e ^* 34 
eauation (3) is riven by frequency, f, samples per second, and converted to a 

equation {i) is given oy. fonnat via ADC ^ ^ vaJues ^ 

35 loaded into a k* point real valued array, called the re- 
xnk) = 0 k = o, k x n (5) ceive time domain array (TDA), in DFT block 31. A 

and DFT is computed of the array to obtain the complex 

valued frequency domain array (FDA) of the sampled 

JWJtt = U A,^xn(-i&n\ \ <6) Af\ in P Ut 011 Kne 3Z ^ k * P° ml DFT ^ COmpleX 

- k) = lUiJmik) r s * s * x/2 " 1 valued receive FDA - Yl(k) ' comainin «' in its first half ' 

' the amplitude and phase modulation information, Au 

„ . . and $\k f of all the harmonics of the transmit signal 

The discrete ttme signal defined by equation (3) is waveform. 



&)$^ k ^^!"ll. V /" e DiSCretC F ° Urier TnmS ' as If the non-linear nature of ADC 29 is ignored, the 
~" "~ demodulation operation can be considered linear. Only 

K out of the k x complex values are retained for decod- 



form (IDFT) of equation (5): 45 

a)(n)-IDFT[Jfi(Jk)]; 0S*./iS* x -l. (7) 



ing. Also, the upper half of the receive FDA is redun- 
To summarize, the 1 baud 15 is generated by taking ix complex conjugates of the values at 

the IDFT of a complex valued array of length k, called 50 Ae image frequencies m the lower half of the array and 
the frequency domain array (FDA). The first half.of the * not retame4 for processmg. Also non-used harmonics 
array is loaded with the amplitudes and phases of the m lowc f half of * e m J ^ 1 ^. the of 
tones to be included in the MFM signal at the corre- noise > °° nXmn MT0 valucs and ^ discarded. When 
spending harmonic numbers. The second half of the ^ do noisc valu «> dlscard 15 equivalent 

array is loaded with the complex conjugate of values in 55 *> filtering out the unused portion of the spectrum at the 
the first half of the array at the image harmonics (the om P ut of demodulation stage. Alternatively, the 
image harmonic of k is k x -k). The resulting data struc- noisc values may be used to estimate the noise level in 
ture is shown for k x = 16 in Table 1. The complex conju- ^ e s y stem . 

gate symmetry about the midpoint insures that the M b y equations (8) and (9), derived in a well- 
IDFT, also an array of Y x points called the time domain 60 known manner from equations (2) and (4), the tones 13 
array (TDA), will be real valued. The values in the m a11 orthogonal; in discrete time on the interval 
TDA, xi(n), defined by equation (7) comprise the dis- O^n^-l, and in continuous time on the interval 
crete time samples for the digitally modulated analog O^u^AT. 
MFM transmit signal waveform. 

TABLE I 65 r Ar <« 
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-continued selected for the transmitted symbol; thus, the symbol 

(9) error will contain only one bit error. Gray encoding 

k *~ 1 ^wn) = * * * ' als0 a* 10 ^ 5 decoding directly into di-bits, with the right 

ii-o " \{A k ) 2 k x x k = ( bit determined by the polarity of the real component 

5 and the left bit by the polarity of the imaginary compo- 

Signals in different bauds are orthogonal, since they nent. The digital data is obtained by successively decod- 

do not overlap in (real) time. Thus, the K tones of the L ing each tone of each baud. 

bauds form a set of LK orthogonal signals so that the Conventional QPSK requires strict phase coherent 
response of a matched filter, or correlator, to any of the regeneration of the sampling frequency to eliminate 
tones (regardless of its modulation) other than the one 10 phase ambiguity. This results in a complex synchroniza- 
to which it is matched is zero. Since the receiver 23 is tion design or distribution of the multitude of carrier 
linear, and a matched filter assures that it is optimal for frequencies to each receiver. MFDQPSK provides a 
demodulation of MFM, the signal-to-noise ratio is maxi- practical solution to the phase uncertainty problem by 
mized in the presence of additive white noise. This is so translating the original symbol set into a second "differ- 
in spite of unknown communication channel or line 15 cnt iaJ" symbol set, or set of states which can be evalu- 
attenuation at each of the frequencies. The channel need atec j to provide a set of modulation values, which is then 
only have constant attenuation over narrow frequency encoded as QPSK. Translated symbols are generated 
bands corresponding to the bandwidths of the individ- bascd on ^ - mpQt sym bol and the previous translated* 
ual tones, Af, and not over the entire bandwidth W. symbo i; F i G S. 4a-td illustrate this translation. As 
Thus the demodulation is still optimal for a widebancj 20 regardlcss of ^ prcvious 8ymbol m ^ of 
signal propagating through a channel with substantial W genera tes a new symbol in the same quadrant as the 
variation in gain across the band. previous symbol. An input of W rotates the new sym- 
Modulation bol -fir/2 radians from the previous, 'IT rotates ir 
w j i r ^ , - . , * r radians, and '10* rotates — 7r/2. In the receiver, decod- 
Modulation of the MFM signal can be uniquely ac- 25 performed by determining the phase difference 
comphshed in a number of ways. The signal x, n) as J twc ^ uccessiv J paks of tone l P 
defined in equation (3), can be modulated in amplitude, "^ w «=" ^^'^ 

frequency In* phase, by translating the da?a into *» ^absence of noise, the phase difference will be 0 

changes in A U , k and <f>u respectively. However, any ±*£ or ff ^ ians ' ^ r ^. Ae P h ^^^°T 

combination of the modulation types is also possible, 30 QPSK constellation +ir/4radiaiis is added to the phase 

such as changing amplitude and phase to produce quad- differences and the quadrature components of the com- 

rature amplitude modulation (QAM). The MFM signal P le * '^presentation of the result are decoded as with 

can be bandpass or baseband and through further multi- W ^ P°J"* of the real part determines the 

plication with a carrier frequency, Xl (n) can be trans- "* ht h ™ d ^ f> d the polarity of the imaginary part 

lated to any frequency band desired. Bandpass signals 35 determines the left hand bit. 

modulated using quadrature phase shift keying (QPSK) MFM SYSTEM 

and differential quadrature phase shift keying (DQPSK) „. „ m 

are described in detail hereinbelow. Referring now to FIGS. 5-8, FIG. 5 is a functional 

Conventional QPSK converts a digital input into four block diagram of an MFM commumcations system 

modulation voltage levels (symbols) to determine the 40 according to the principles of the present invention, 

phase of the transmitter output. To minimize the proba- The MFM communications system of the present inven- 

bility of symbol error the phases are spaced at multiples tion comprises a transmitter section 51 wherein digital 

of ir/2. A plot of the complex envelope of one tone of **** m P ut on 62 » modulated in a desired frequency 

x,(n) is shown in FIG. 3. The angle <f> can take on the band to provide an MFM analog transmit signal which 

values of ±7r/4 and ±3?r /4. 45 * s coupled to a communications channel 55 for propaga- 

Encoding MFM with QPSK is accomplished by tion to a receiver section 53. The receiver section 53 

loading the complex frequency domain array with the demodulates the received analog signal including any 

appropriate phase information. For example, if the digi- n °ise introduced in the communication channel 55 to 

tal input is *0ll0 . . . \ The first di-bit or symbol, '01', provide the original source data output on line 66 in 

would be loaded into the frequency bin, k, as RepC(k)]- 50 digital form. The input data on line 62 may be any infor- 

= - Ajfckjc/2(2)i and Im[X(k)]=: - Re[X(k)] where A* is mation or message which may be represented in a digi- 

the amplitude of tone k. Similarly, the second symbol tal format. While the modulation and transmission pro- 

'10* would be loaded into the second bin, k+ 1, as cess is preferably in real time, input information can be 

Re[X(k+ 1)]= A* + ikx/2(2)* and Im[X(k-h 1)]=- processed and stored as a time domain array (TDA) in 

— Re[X(k+ 1)]. This continues until all the tones in the 55 the system memory for transmission at a later time. All 

baud are filled or there are no more symbols. As de- modulation/demodulation and encoding/decoding 

scribed above, the discrete time domain signal is pro- functions are performed entirely internally by industry 

duced by taking the IDFT of the complex frequency standard microcomputers, such as personal computers 

domain array. Each successive baud is encoded simi- (an IBM AT Personal Computer including a floating 

larly. 60 point array processor is suitable for this application), 

Decoding QPSK directly into bits is accomplished under the control of special purpose transmit and re- 

easily as follows. Assuming a coherent receiver, decod- ceive computer programs thereby eliminating the need 

ing requires evaluation of the polarity of the real and for special purpose modems to translate between the 

imaginary components of the K complex values stored digital and analog domains. The transmit microcom- 

in the receive frequency domain array. As shown in 65 puter 57 is coupled to the communication channel 55 via 

FIG. 3, the symbol mapping uses Gray encoding. This its I/O ports 63 and DAC 65. The receive microcom- 

reduces the probability of bit error as errors caused by puter 59 is coupled to the communications channel 55 

AWGN are likely to occur when the adjacent symbol is via ADC 67 and its internal I/O ports 71. Matched filter 
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69 provides baud synchronization within each packet. ceiver system with a real time and a non-real time DFT 

The communications channel 55 may be any available processing capability, 

communications link linking the transmit 57 and the TABLE II 
receive 59 microcomputers together. The medium for 
transmitting the data may be wire or optical fiber, or it 5 
may be radio frequency such as HF, microwave or 
satellite. 

FIG. 6 is a functional block diagram of the transmit- 
ter section 51. As described above, the input data is 
processed and may be transmitted immediately, i.e., in 10 
real time, or stored for transmission at a later time. The p or tn e MFM system of the present invention, carrier 
digital input to the encoder 511 is a serial string of a frequency synchronization is not required because of 
binary digits. Processing by the encoder 511 is per- tnc differential encoding. In the preferred embodiment, 
formed on a baud by baud basis until the end of the data the MFM signal is implemented in the form of a trans- 
file or the maximum packet length has been obtained. 15 mission packet. To acquire the packet, the receiver 
The encoder 511 encodes the input data words or sym- section 53 must know when to start sampling. This is 
bols into complex valued symbols which are stored in a accomplished by packet synchronization. Typically, 
frequency domain array (FDA), a file defined in the unique words are inserted to mark the start of each 
system memory (not shown). The value of the input packet or frame. In the preferred embodiment, the 
symbols is a function of the type of modulation. For unique word is termed a "synchronization baud", and is 
example, QPSK symbols correspond to two-bit, or di- added at the beginning of each packet The synchroni- 
bit, words and are encoded as discussed hereinabove. A zation baud is generated similarly to the other bauds in 
discrete signal is produced by computing the IDFT 513 the packet, except that the tones and phases are prede- 
of the complex valued FDA array and is loaded into the 2J termined. Acquisition of data by the receiver section 53 
packet storage area 515 in a time domain array (TDA). commences after successful detection of the synchroni- 
The controller 519 determines the parameters of the zation baud. 

packet based on the modulation type, baud rate and FIG. 6 is a detailed block diagram of the data acquisi- 

'message size. The controller sequences the input data tion board 531 for the receive section 53. Frame syn- 

through the transmitter section 51 one baud at a time ^ chronization of DQPSK MFM in the present invention 

until the entire message has been processed. Once pro- is achieved utilizing a polarity-only matched filter or 

cessed, the entire stored digital signal, xi(n), is trans- polarity coincidence correlator. The data acquisition 

ferred from the TDA at a selected clock rate, f x samples board 531 comprises ADC 67 which couples a 12-bit 

per second, through DAC 517. Depending on the com- digital data signal on line 541 to AND gate 547 and also 

munication channel 55 frequency response, filtering of 35 couples a single bit signal on line 542 to the matched 

the analog output signal may be desired to remove high filter 69 (as shown in FIG. 5). The matched filter 69 

frequency amplitude discontinuities introduced by provides sampling triggers to the AND gate 547 syn- 

DAC 517. chronized with respect to the time of arrival of the input 

FIG. 7 is a functional block diagram of the receiver data packet. Using the polarity of the analog input, 
section 53. The data acquisition board 531 samples the 40 coupled as a single bit on line 542, it functions as a 

input analog signal at the clock rate f x samples per sec- matched filter to the last half of the 256 point synchroni- 

ond, and converts it to a digital format via ADC 67. zation baud. The hard limiter 544 comprises a fast, high 

Filtering of the analog input signal may be required to precision, high gain operational amplifier utilized to 

bandlimit input noise and to reject out of band interfer- obtain the polarity information in the analog input sig- 
ence from frequencies greater than fx/2. The digital 45 nal. During the positive and negative portions of the 

data sequence thus obtained is stored as a time domain one-bit digital input on line 542, the output of the hard 

packet or array (TDA) in the receive microcomputer 59 limiter 544 is positive and zero, respectively. This uni- 

data storage 533. The stored real value data sequence is polar signal is synchronized to the receiver section 53 

accepted from the receive TDA one baud at a time to clock frequency, f x , as it is clocked through a 128 point 
perform a k x point DFT 535. The first half of the result- 50 polarity-only matched filter 543, To minimize bit insta- 

ing complex values are decoded in decode block 537 to bility the hard limiter 544 slew rate should be as large as 

obtain the amplitude and phase modulation information. possible, (an LM 318 general purpose operational ampli- 

In a manner similar to the transmit section 51, the re- fier having a slew rate of 70 volts per microsecond is 

ceiver section controller 539 sequences the data suitable for this purpose.) 

through the receive section 53. 55 In the preferred embodiment, the polarity-only 

As a function of the processing speed and capacity of matched filter 543 described above is realized by a com- 

the receive microcomputer 59, some delay may be in- mercially available CMOS programmable digital corre- 

troduced to the receive section from input to output. lator comprising four, 32-bit correlator modules in se- 

For real time processing, the receive microcomputer 59 ries to form a 1 X 128-bit fixed channel configuration. A 
must have the capability to complete the kx point DFT 60 single bit is clocked in from ADC 67 on line 542 to be 

535 in the same or less time than it takes to fill up one compared against a 128-bit preloaded reference pattern, 

baud length buffer in the receive TDA. In a microcom- A match occurs when the synchronization baud for 

puter utilizing parallel TDP buffers, data processing can each packet in the received signal is clocked through 

alternate between buffers and data flow will be continu- ADC 67 (a TMC2221 TTL compatible CMOS correla- 
ous through the system at the sampling rate. The stor- 65 tor manufactured by TRW LSI Products Inc. is suitable 

age requirements and processing time per packet for a for this purpose.) A synchronization trigger generated 

clock rate of f x =61,440 Hz, L=30 bauds in a packet when a match occurs enables D flip-flop 548 to couple 

and k x =1024 samples are given in Table II for a re- sampling triggers to AND gate 547 for data acquistion 
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at the beginning of each data packet as required for the 
demodulation process. 

In the preferred embodiment, the MFM communica- 
tions system of the present invention may be imple- 
mented in a variety of ways utilizing various software 
and hardware configurations. To provide maximum 
I/O channel throughput, the host computer direct 
memory access (DMA) is utilized. Any desired packet 
construction may be used and various baud rates are 
available to provide greater flexibility in adapting the 
packet to available channel parameters. For example, in 
MFDQPSK modulation, where the information is rep- 
resented by the phase difference between adjacent 
tones, channel tone distortion affects a shorter baud 
more than a longer baud due to the greater tone spac- 
ing. A longer baud is preferred when differentia] encod- 
ing is between adjacent tones. A short duration baud is 
desirable with differential coding between the same 
tones on adjacent bauds. 

Referring now to FIG. 9, a two-stage, data driven 20 
echo canceller (EC) 70 according to the present inven- 
tion makes use of a unique and important feature of 
MFM, namely that data is available during the en- 
coding/modulation process and during the signal trans- 
mitted demodulation/decoding process both as digital 25 
time domain sequences and as frequency domain arrays. 
The first stage of the EC 70 uses a finite impulse re- 
sponse (FIR) linear filter h(n) 79 driven by the transmit 
time domain digital signal sequence x(n) output by the 
IDFT 75 to generate a first estimate of the echo signal. 
This first echo estimate is coupled via DAC 81 to sum- 
mer 86 at the channel input 78 and subtracted from y(t), 
the received analog signal prior to its analog to digital 
conversion. The pre-con version cancellation of the 
received echo reduces the dynamic range requirement 
for the receive ADC 85 to a maximum of 14 bits. The 
first stage utilizes the same x(n) data which is transmit- 
ted on the channel 76 via DAC 77. 

The second stage of the EC 70 comprises multiplier 
85 and summer 87 and operates in the digital frequency 40 
domain using the transmit frequency domain array X(k) 
72 as the data driver. The complex values in the trans- 
mit frequency domain array are multiplied by complex 
coefficients C(k) at multiplier 85 to generate an estimate 
of the residual echo that remains after the first stage of 45 
cancellation. This refined estimate is subtracted from 
R(k), the values in the receive frequency domain array 
83. The second stage of EC 70 cancels residual near-end 
and far-end echo whereas the first stage of the EC is 
aimed primarily at reducing the near-end echo. 50 

The coefficients h(n) and C(k) for both stages of the 
EC are established using a training packet prior to the 
beginning of data exchange. The algorithms for com- 
puting the coefficients are described below. The two 
nodes are trained sequentially. That is, the call initiating 55 
modem transmits a training packet that is utilized to 
adapt its echo canceller with the answering modem 
connected, but not transmitting. Then, the answering 
modem trains its EC with the calling modem connected 
but not transmitting. Any other echo cancellers or echo 60 
suppressors in the link are disabled. After both ECs are 
trained, full-duplex data exchange proceeds using the 
aforementioned differentially encoded MFM signals. 

Referring now also to FIG. 10, algorithms are de- 



tl. 30 



35 



ond stage of the EC is trained using the remainder of the 
bauds in the training packet with the first stage EC in 
operation. Thus the training packet must contain at least 
two bauds but may contain more to improve the esti- 
mate of the residual echo. In one preferred embodiment, 
the training packet contains 8 bauds. Two embodiments 
for computing the N filter coefficients are given by 



Iterative Algorithm 

A-+ l © « *"(0 + k~ l (/{n)x{n - 0); 



00) 



. 0, 1 . . . N - 1; n - 0, 1 . . . ft* - 1; 



and 

Smoothed Algorithm 

n=0 n=0 



01) 



im 0, 1, ... N - 1. 

In all cases, N, the number of coefficients in the first 
stage echo canceller, must be very much less than k x , 
the number of sample points in a baud. In one preferred 
embodiment, N= 1 and kjc=1024. However, all FIR 
filters used in conjunction with MFM digital data trans- 
mission echo cancellation. Both equations (10) and (11) 
are LMS (least mean square) type algorithms and are 
effective for echo cancellation; equation (11) is more 
computationally intensive than equation (10) but in 
general results in slightly more echo reduction. 

Referring now to FIG. 11 and the second stage of 
EC, two embodiments for computing the K complex 
weighting coefficients are given as 



Iterative Algorithm 

d+»<A) = d(k) + L-*[RKk) xmi 

I » 1, 2, ...£.-*» ft], ... 
and 



02) 



Smoothed Algorithm 



ftjr-1 ftx-1 
■ 2 *i(ft)Wft)/ 2 
n=0 rsO 



03) 



Again, both algorithms are effective; however, in order 
to keep the training packet as short as possible, L, the 
number of bauds may be minimized. When there are 
very few bauds, the iterative algorithm will not be as 
desireable as the smoothed algorithm unless the initial 
guess in the iterative algorithm is very nearly correct. 
In one preferred embodiment, equation (10) is used with 
equation (13) and N=l and L=l; while in another 
preferred embodiment, equation (11) is used with equa- 
tion (13) and N=2 and L=2. 

The echo canceller of the present invention provides 
up to 75 db of echo rejection for MFM communica- 
tions. This makes possible full-duplex signalling over 
dial-up two-wire telephone circuits at bandwidth effi- 
ciencies up to 8 bits per Hz. The EC may be imple- 



scribed for training each stage of the EC. The first stage 65 mented using commercially available 14 bit ADCs and 



of the EC is trained using the first baud of the training 
packet. As shown in FIG. 1, a baud 15 is of length AT 
seconds and contains k x digital signal samples. The sec- 



DACs, for example, the Texas Instruments TIC32044. 
Two stages of LMS type echo cancellation are used to 
obtain this level of performance; stage one operates in 



07/29/2004, east version: 1.4.1 



17 



5,166,924 



18 



the time domain while stage two operates in the fre- 
quency domain. Both stages are driven by the transmit 
data. 

MFDQPSK Modulation 

Referring now to FIGS. 12, 13 and 4a-4d, FIGS. 12 
and 13 are a flow chart illustrating multi-frequency 
differential quadrature phase shift keyed (MFDQPSK) 
modulation and demodulation in which 2 bits of data are 
encoded into each tone. The resultant MFM signal 10 
transmits 2 bits per Hz of the channel bandwidth occu- 
pied. 

The encoding and generation of the MFDQPSK 
signal is illustrated in FIG. 12 in which all blocks 
91-101 may be resident in a computer software program IS 
executed by the transmit microcomputer 57. Alterna- 
tively, some or all of the blocks may be implemented in 
hardware or firmware, such as in a DSP chip. Input 
data to be transmitted are placed in a Data File 1 (91) in 
a 2-bit or di-bit format. The differential encoding step 93 20 
differentially encodes the phase transferring the en- 
coded di-bits to Data File 2 (95). Data File 2 is created 
with an additional initial state set to So- The first di-bit 
from Data File 1 is then encoded into the second state of 
Data File 2 according to the state transition diagram 25 
shown in FIG. 4c. The second di-bit from Data File 1 is 
encoded into the third state of Data File 2 in accordance 
with the state transition diagrams of FIGS. 4a-4d where 
the beginning state is the previous state of Data File 2. 
This process is continued until all K di-bits from Data 30 
File 1 have been differentially encoded into K+ 1 states 
of Data File 2. Each of the K-f 1 digital states provide 
a modulation value for modulating one of the set of 
available tones. Examples of Data File 1 and Data File 
2 are shown in Table III. 35 

Encoding Step 2 (97) provides amplitude and phase 
encoding of the K + 1 states from Data File 2 into the 
available tones or carrier frequencies and creates a com- 
plex-valued frequency domain array (FDA) file 99 con- 
sisting of k, complex entries. The entries in both Data 40 
File 1 and Data File 2 and on lines 92 and 94 comprise 
n-bit words (n=2 for this embodiment). K+ 1 of the 
complex entries in this file somewhere between the first 
entry, entry 0, and entry kjc/2 — 1 will be encoded ac- 
cording to the states of Data File 2 as follows: Assume 45 
entry number kt is the first non-zero entry in the FDA 
File 99. The first state in Data File 2 is always S 0 so 
entry ki will be loaded with a real part of 1. and an 
imaginary part of 1. Subsequent entries in the FDA 99 
are loaded with complex values according to the subse- 50 
quent data stored in Data File 2 according to the fol- 
lowing rule: If the least significant bit of the state index 
expressed as a binary number is logical zero, the real 
part r is set to 1.0; if the least significant bit is logical 
one, the real part r is set to — 1. If the most significant 55 
bit of the state index expressed as a binary number is 
logical zero, the imaginary part i of the entry is set to 1.0 
and if it is logical one the imaginary part i is set to — 1. 
This procedure continues until entry k2=ki-fK is 
reached and all data in Data File 2 have been encoded 60 
into complex valued digital signal values. The remain- 
ing entries in the first half of the FDA 99 are set to 0. 
The second half of the entries in FDA are encoded 
according to the complex, conjugate image rule. This 
rule is the following: The image of entry k in the first 65 
half of the FDA 99 is entry k-c-k. The real part of entry 
k*— k is set to the same value as the real part of entry k. 
The imaginary part of entry k x — k is set to the negative 



of the imaginary part of entry k. The procedure de- 
scribed above for the Differential Encoding Step 93 and 
the Amplitude and Phase encoding step 97 is illustrated 
in Table III for a Data File 1 consisting of 4 di-bits and 
in Table IV an FDA 99 consisting of 16 complex 
entries. In this illustration, k\ has been chosen to be 
entry 2. 

TABLE III 



Data 

File 1 



Data 

File 2 



01 
10 
11 
00 



So 
S| 
So 
S 3 
S 3 



TABLE IV 




Frequency Domain 






Array 




k 


i 


r 


0 


0.0 


0.0 


1 


0.0 


0.0 


2 


1.0 


1.0 


3 


1.0 


-1.0 


4 


1.0 


1.0 


5 


-1.0 


-1.0 


6 


-1.0 


-1.0 


7 


0.0 


0.0 


8 


0.0 


0.0 


9 


0.0 


0.0 


10 


1.0 


-1.0 


11 


1.0 


-1.0 


12 


-1.0 


1.0 


13 


-1.0 


-1.0 


14 


-1.0 


1.0 


15 


0.0 


0.0 



The FDA 99 is transformed to the time domain array 
(TDA) 101 by means of a k x point Inverse Discrete 
Fourier Transformation (IDFT) 100. Because of the 
complex conjugate image symmetry introduced in the 
FDA 99 the result of the IDFT 100 is real-valued. Thus 
the TDA 101 consists of k* real-valued entires. Each 
real value is stored as an integer-valued byte. These 
bytes are sent via an I/O port 63 of the transmitting 
computer 57 to a DAC 65 which is clocked by a clock 
64 operating at f x clock pulses per second. An external 
clock may be provided. The clock rate f Xt may be any 
value from 0 to an upper limit set by the host computer's 
57 ability to transfer bytes from the TDA 101 to the I/O 
port 63. In a particular embodiment f*= 61440. samples 
per second, the host computer comprises an IBM 
PC/AT and the I/O port utilized is the direct memory 
access (DMA) channel. The DAC 65 converts the byte 
representation of the digital signal samples to voltage 
levels thus providing an analog transmit signal which is 
coupled to the communications channel 55 connecting 
the data source and destination. The 1c* TDA bytes 
representing the discrete time samples of one baud are 
clocked out of the DAC 65 in an interval AT=kx/f* 
seconds. Therefore the clock rate of the system, i xt is k* 
times the baud rate, 1/AT. 

The analog transmit signal contains K+l tones 
spaced at intervals of 1/AT beginning at frequency 
ki/AT and ending with frequency ki/AT. The maxi- 
mum value for k2 is (kx/2— 1) which is less than f*/2 as 
required by the Nyquist Sampling Theorem. The 
MFDQPSK modulation system operates at a baud rate 
of 1/AT with a bit rate of 2K/AT. In another embodi- 
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ment, the system implemented selectively operates at 
baud rates of 15, 30, 60, 120, and 240 bauds per second 
with a constant bit rate of 8192 bits per second in a 4000 
Hz bandpass channel. 

The signal flow diagram for receiving and decoding 5 
the MFDQPSK signal is shown in FIG. 13. The re- 
ceived analog signal is synchronized and digitized by 
data acquisition board 531 (as shown in FIGS. 7 and 8) 
which couples digitized, real-valued data to the receive 
time domain array or packet (TDA) 103. Process blocks 10 
103-115 are internal to the receiving computer and are 
implemented as a computer software program. Alterna- 
tively, some or all of process blocks may be imple- 
mented in hardware or firmware. Receive ADC 67 
must be clocked at the same rate, f Xl as the DAC 65 in 15 
transmit section 51 with an accuracy of plus or minus 
fx/(16 ki) Hz. Baud synchronization timing accuracy 
must be established to within plus or minus four clock 
cycles. In the preferred embodiment this is established 
without direct timing information from the transmit 20 
section 51 by preceding each packet of information 
bauds with a known synchronization baud and correlat- 
ing it with a stored reference in matched filter 69 as 
described hereinabove. 

The data acquisition board 531 provides digitized 25 
received signal samples to a receive computer I/O port 
71. The host receive computer 59 stores these samples 
in the receive TDA 103 having k x times L entries in a 
real-valued array. The DFT 105 and Decoding Step 1 
109 transform k x entries at a time, that is one baud at a 30 
time, from the receive TDA 103 to the complex valued 
receive digital signal array (DSA) 111 according to the 
following procedures: The k* point DFT of the first k* 
real valued points in the TDA 103 is obtained and only 
those (complex valued) coefficients between entry ki 35 
and k2 are retained and placed in a temporary digital 
signal array 107. All other components of the transfor- 
mation are discarded. The entry corresponding to ki -f 1 
is multiplied by the complex conjugate of the entry 
corresponding to kt to create a complex number with a 40 
phase equal to the phase difference between tones ki + 1 
and ki. This complex number then is multiplied by the 
complex number exp(j7r/4), to add ir/4 radians to the 
phase difference, and the result is stored as the first 
entry in the receive DSA 111. The same operation is 45 
performed on entries k"i +2 and ki + 1 and so forth until 
the final operation is performed on entries k2 and k2 — 1 
resulting in K complex entries in the receive DSA 111. 
The K complex entries are decoded in Decoding Step 2 
(113) to create Data File 3 (115) as follows: if the real 50 
part of the complex entry is positive, the least significant 
bit of the di-bit corresponding to that entry in the DSA 
111 is made logical zero; if the least significant bit is 
negative, it is made logical one. The most significant bit 
of the di-bit is decoded from the imaginary part of the 55 
corresponding complex entry in the same manner. The 
final result is that the first K di-bits of Data File 1 (91) 
at the source (i.e., input data at the transmit section 51) 
have been transferred to Data File 3 (115) at the destina- 
tion (i.e., receive section 53). This procedure is repeated 60 
on the remaining entries in the real valued data stored in 
the receive TDA 103, k* values at a time, until all L 
bauds in the packet have been decoded. In the preferred 
embodiment, with k* of 1024, a baud rate of 60 bauds 
per second, K of 64 and L of 30 bauds, 3840 bits of 65 
information are transmitted in a packet having an inter- 
val of i second in a 3960 Hz bandwidth. The band may 
be placed anywhere in the frequency spectrum from 60 
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Hz to 30,660 Hz simply by choosing different values for 
ki. 

MFD16QAM Modulation 

Referring now also to FIGS. 14c, 14* and 15, in the 
preferred embodiment, multi-frequency differential 16- 
quadrature amplitude modulation (MFD16-QAM) dif- 
fers from the MFDQPSK modulation described above 
in the following ways. The digital input data is stored in 
Data File 1 (91) in 4-bit words referred to as "nibbles" 
and are encoded into each tone using the 16 point state 
diagram or constellation shown in FIG. 14a. Thus, in 
MFD16-QAM, data is transmitted at 4 bits per Hz of 
allocated channel bandwidth. 

The state diagram for the differential encoding from 
frequency to frequency is given in Table V for each of 
the 16 values for an input nibble. The diagram defines 
the state for tone k+ 1 for the conditions that tone k is 
in an even numbered state, and that tone k is in an odd 
numbered state. Index n takes on all values from 0 
through 7 and the indicated index n addition is per- 
formed modulo 16. 

TABLE V 



State for tone k + 1 
Given that tone k is in state: 



Input Nibble 


$2» (Even) 


S2* + i(Odd) 


0000 


S2A+1 


s ta 


0001 






0010 




S2ji+2 


0011 




$2n+3 


0110 


$2*+5 


S2n+4 


0111 


&2n+4 


S2ji+3 


0100 


S2ii+7 


S2A+6 


0101 




S2n+7 


1100 


S211+9 


S2n+8 


1101 


S2n+8 


S2b+9 


1110 


S2R+11 


S2n+10 


mi 


S2n+10 


S2n+ll 


1010 


S2b+13 


S2n+12 


1011 


S21H-U 


Sin +13 


1000 


S2n+lS 


S2o+14 


1001 


Sjn+14 


S2n+15 



The first 3 bits i.e., the 3 most significant bits, of each 
input nibble represent phase and are differentially en- 
coded as a Gray code so that differentially decoded ; 
phase symbols that are adjacent in phase only differ by 
1 bit as shown in Table V. The fourth bit Le., the least 
significant bit, of each input nibble is an amplitude bit 
and is encoded such that the magnitude of tone k+ 1 is 
changed to the alternate amplitude if the amplitude bit is 
a logical zero but remains unchanged if it is a logical 
one. Table VI is an example of a 4-nibble message stored 
in Data File 1 and the differentially encoded data in 
Data File 2 (95). Table VII illustrates the corresponding 
FDA file 99 after the amplitude and phase encoding 
process 97. 



TABLE VI 



Data File 1 


Data File 2 




So 


0111 


S4 


1100 




1010 


s 8 


1001 


S6 
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TABLE VII ^ dcs for P acket size » m ^ generates a synchronization 

baud. 



U 


u.u 


u.u 


1 


0.0 


0.0 


2 


.2296 


.5543 


3 


.5543 


-.2296 


4 


-.9239 


.3827 


5 


-.2296 


-.5543 


6 


.2296 


-.5543 


7 


0.0 


0.0 


8 


0.0 


0.0 


9 


0.0 


0.0 


10 


-.2296 


-.5543 


11 


.2296 


-.5543 


12 


.9239 


.3827 


13 


-.5543 


-.2296 


14 


-.2296 


.5543 


15 


0.0 


0.0 



frequency domain array EncodeQPSK 163 begins by displaying a four symbol 

QPSK constellation which is used as a reference in 
5 selecting the symbols over the band. The symbols for 
the K tones in the baud can be selected in a variety of 
ways: symbols for all tones may be randomly selected 
from a random generator; a symbol for each tone may 
be selected by the user; and individual tones may be 
10 removed from the band. The capability to remove tones 
from the band allows construction of a baud with an 
arbitrary number of tones within the band determined 
by ki and ki. Selected symbols for the band are loaded 
into the complex frequency domain array with their 
15 complex conjugate image frequency. EncodeQAM 165 
is functionally identical to EncodeQPSK 163, except 
the symbols conform to a 16-QAM constellation (as 
shown in FIG. 14a). 

_ , , . r . wTTT^ir ^ a w tt\d To obtain the real discrete time-domain array of the 

Demodulation of the ^16-QAM receive TOP ^ encoded baud, ComplexFFT 167 computes the inverse 
proceeds similarly as demodulation of the MFDQPSK ^ ComplexFFT 167 consumes the majority of the 
to the poinl t of ^creation of the temporary digital signal ^ m ^ ro ^ due to thc 

array 107 . A 12-bit ADC is utuaed m arithmetic operations required, thus restricting the 

tion board 531 to reduce quantization noise. To decode ovcral] th h t ofthc system . For rcal ^ proccss . 
the 3-bit phase symbols, each complex entry in the tern- 25 the ^ orithm ^ computed by a floating point 
porary digital signal array 107 is multiplied by the com- ^ processor 

plex conjugate of the previous complex entry and the Each value m' the toe domain array is represented as 
product is multiplied by the complex number expO^r/8). a rea] data ^ occupying six bytes of memory. Scale 
In the absence of noise or other errors, the resulting Data x69 converts these values down to a one byte 
complex products will be the 24 possible points in the 30 fonnat acceptable t0 the DA C and places them into a 
phase diagram or constellation shown m FIG. 14b. The ^ acket storage ouffer EncodeQPSK 163, ComplexFFT 
largest amplitude points occur when adjacent tones 167 ^ 169 m cxe cuted for each baud, 

both have the large amplitude. The smallest amplitudes m ^ ^ L baud have processed. To transmit the 
occur when adjacent tones both have the smaller ampli- packet out of ^ compter, DMAINIT 171 transfers 
tude. The intermediate amplitudes occur when the adja- 35 ^p^ at Ac clock f reqU ency, f x samples per second 
cent tones have different amplitudes, i.e., the amplitude over Ae D MA channel to the DAC. The source code 
bit is a logical zero. listing for the TRANSMIT program is given in Appen- 

The eight possible phase symbols are decoded in dix A attached hereto, 
accordance with the decision tree shown FIG. 15. The Referring now to FIG. 17, a procedural flow diagram 
phase decoding rules shown are amplitude independent 40 of another transmit computer software program named 
and are the optimum decision rules for decoding the XMITMES as implemented in another preferred em- 
phase in the presence of additive noise. bodiment is shown. The XMITMES program generates 

The amplitude bit is decoded by comparing the mag- m mfM data signal which encodes and transmits an 
nitudes of the complex values in adjacent positions of ASCII file using DQPSK. A synchronization baud is 
the temporary digital signal array 107. In the preferred 45 attached at the beginning of each packet. The 
embodiment, the large and small magnitude tones are XMITMES program has a less complex structure than 
created in a ratio of 5 to 3. If the received magnitudes of the TRANSMIT program described above because all 
adjacent tones differ by more than i the magnitude of tones are placed in a fixed band and are encoded utiliz- 
the sum, the amplitude bit is decoded as a logical zero, m g DQPSK. 

otherwise it is decoded as a logical one. 50 The synchronization baud for XMITMES is a prede- 

Upon completion of the decision process in FIG. 15, termined sequence generated by SyncBaud 173. This 
K nibbles or K 4-bit words have been transferred from baud is constructed as are all other baud, except k, is 
Data File 1 (91) at the transmit computer 57 to Data fixed at 256, and tones 68 to 83, which are the tones in 
File 3 (115) of the receive computer 59. This demodula- the fixed band, are encoded with the same random sym- 
tion and decoding process is repeated for each of the L 55 bol pattern regardless of the packet construction or 
bauds in the receive TDA 103 until 4KL bits (KL/2 input message 176. This synchronization sequence oc- 
bytes) have been transferred from Data File 1 at the cupiesthe first 256 values in the packet buffer and there- 
transmitter 51 to Data File 3 of the receiver 53. In the fore is the first to be clocked out of the transmit corn- 
preferred embodiment with k x of 2048, K of 128, L of 15 puter 57. 

and baud rate of 30 bauds per second, 960 bytes are 60 Since the input message 176 can be of arbitrary size, 
transferred in J second in a 3960 Hz bandwidth at TailorPacket 174 determines the maximum number of 
15,360 bits/sec. baud required for encoding. This is determined by di- 

Referring to FIG. 16, a procedural flow diagram of a viding the number of characters in the message by the 
computer software program named TRANSMIT as number of characters that can be encoded, 
implemented in the preferred embodiment is shown. 65 DiflEncode 177 encodes the message file 176 into the 
The TRANSMIT program generates the MFM data complex frequency domain-array. It reads one charac- 
signal and provides encoding of data files utilizing 16- ter at a time; then breaks the 8-bit character into four 
QAM, QPSK and DQPSK methods. The program pro- 2-bit words. The words are DQPSK encoded and 
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stored in the frequency domain array. Once encoded, 
processing and signal output by ComplexFFT 179, 
ScaleData 181, and DMAINIT 183 proceed similarly as 
in TRANSMIT program. The source code listing for 
the XMITMES program is given in Appendix B at- 
tached hereto. 

Referring now to FIG. 18, a procedural flow diagram 
of a computer software program named RECEIVE as 
implemented in the preferred embodiment is shown. 



ization, conversions take place on the positive transition 
of every trigger until the specified number of samples 
have been acquired and transferred to memory. DhT- 
Decode 191 demodulates the encoded symbols by dif- 
ferentially decoding the complex frequency array trans- 
formed from the sampled data by ComplexFFT 189. 
AcquireData 187 samples and stores a memory segment 
of data regardless of the packet si ze, L k x - However, 
once stored, ConvertData, ComplexFFT 189, and Diff- 



The RECEIVE program demodulates the encoded 10 Decode 191 process only L baud of the data. Diff- 



MFM transmit signal generated by the TRANSMIT 
program. A user interactive procedure, PacketSetUp 
185 tailors the processing to the expected receive 
packet, using the inputs k x and L. AcquireData 187 



Decode 191 combines four differentially decoded sym- 
bols into one byte, representing the ordinal number of 
an ASCII character. To reconstruct the message, the 
characters are transferred to text file MESSAGE.DAT 



allocates memory to store the sampled values trans- IS until processing is complete. When desired, ShowMes- 

ferred from the board using the receive computer DMA sage 193 displays the recovered message. The source 

controller. Direct Memory Access is the only data code listing for the RECEIVE program is given in 

transfer mode capable of transferring data to memory at Appendix C attached hereto. 

the required rate without disruption by other interrupt While the present invention has been particularly 

processes in the computer! Other important initializa- 20 shown and described with respect to certain preferred 

tion parameters are triggering source and the number of embodiments thereof, it should be readily apparent to 

samples to be collected. The ADC may be triggered those of ordinary skill in the art that various changes 

from two sources, a programmable interval timer or an and modifications in form and details may be made 

external trigger source (not shown). The programmable without departing from the spirit and scope of the in- 

interval timer divides either a 1 MHz or 10 MHz clock 25 vention as set forth in the appended claims, 
to derive the sampling rate of the trigger. After initial- 



. APPENDIX A TRANSMIT 

program TRANSMIT; 

(transmits variable length packots up to 61440. Tho length is 
determined by BCSTARRAY and can bo tako to a lull. segment of memory. 
No sync baud is uaod. All external flics nro aanumod to bo in Hie r.nma 
directory*) 



type 

TNvector * array[ 0. . 4095] of real; 

TNvectorPtr = ATNvector; 

BCSTARRAY » array{ -28673. • 32767J of byte; 



(*stores output samples*) 



var 

kx, (*baud size*) 

kl, (* lowest tone in band*) 

k2, (*highest tone in band*) 

NUHBAUDS, (*numbar of bauds in packet*) 

BAUD COUNT, (*curr ant baud being processed*) 

BYTEC0UNT, (*number of bytes to be transferred by DMA*) 

I : integer; 

XREAL, (*complex freq/time domain arrays*) 

XIMAG : TNvectorPtr; 

QAM, (*encoding scheme*) 

INVERSE -.boolean; (*direction of FFT*) 

ERROR :byte; (*status of ComplexFFT*) 

BCST : BCSTARRAY; (*packet storage buffer*) 

T3IEFILE : f ile of byte;(*time sequence output file*) 

ANSWER :char; (*input variable*) 

(*$I FFT87B2. INC*) 
(*$I C0MPFFT. INC*) 

(* - *) 

procedure SelectBaud; 

(*SelectBaud establishes QPSK or QAM encoding kx, kl, k2 9 and the number 
of baud in the packet*) 
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var 

ANSWER : integer; 
MODANSWER : char; 

begin 

vindow( 20,10,80, 20); 

(''Select encoding scheme*) 
repeat 
clrscr; 

writcln( 'Select modulation desired 1 ); 
vriteln; 

writelnC A --- QPSK 1 )-, • 

writelnC ' B — 16 QAM'); 

readlnC MODANSWER); 
until MODANSWER in I • n 1 . ' A ' , ' b 1 ( ' D 1 ] ; 
if MODANSWER in [ V/b'] then 

QAM :■ true; 

(*Select buad size*) 
kx:=0; 
repeat 
clrscr* 

if kx < 0 then writelnC'TRY AGAIN'); 

writelnC 'What is the length of the bauds (kx)?'); 

writelnC 1 i. e. 256, 512, 1024, 2048, 4096* ); 

readlnC ANSWER); 

case ANSWER of 

256: kx:«256; 

512: kx:«512; 

1024: kx:=1024; 

2046: kx:=2048; 

4096: kx:=4096; 
end; C*case ANSWER*) 
if kx = 0 then kx := -1; 
until kx > 0; 

(*Set bandliraits*) 
case kx of 
256: begin 

kl:=68; k2:=83; 
end; 
512: begin 

kl:-135; k2:=166; 
end; 
1024: begin 

kl:«269; k2:=332; 
end; 
2048: begin 

kl:=537; k2:=664; 
end; 
4096: begin 

ki:=1073; k2:-1328; 
end; 

end; C*case kx*) 

(*Select number of baud*) 

writeC 1 How many bauds do your desire to transmit? '); 

reodlnCNUHBAUDS); 

windowC^l^O^S); 
end; (*SelectBaud*) 
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(* : -*) 

procedure DisplayQAM; 

(*DisplayQAM shows the 16-QAM constellation .in the upper portion of the 
screen. A window is set at the bottom of the screen for further 
interaction with the user.*) 



begin 
clrscr; 
gotoXYCll.l); 
vrite(chr(218)); 
for I:» .12 to 64 do 

write(chr(196)); 
write(chr(191)); 
for I:= 2 to 19 do 

begin 

GotoXY(ll.,I); 
GotoXY(65,I); 

end; 
GotoXYC 11,20); 
write(chr(192)); 
for I : = 12 to 64 do 

write(chr(196)); 
write(chr(217)); 

window( 12,2,63,19); 
write('This program 
writelnC signal. ' ); 

writeln('The vector space is shown for one frequency. 1 ); 



write(chr(179)); 
write(chr(]79)); 



encodes a 16 *QAM multifrequency 1 ); 



writelnC 1 
writeln( 1 
writelnC 1 
writelnC 1 
writelnC 1 
write ln( 1 
writelnC 1 
writeln( 1 
writelnC 1 
writelnC 1 
writelnC ' 
writelnC 1 
writelnC' 
writelnC ' 



3 
7 



12 
16 



11 
15 



,chr(179), 
,chr(179), 
,chr(179) t 
,chr(179), 
,chr(179), 
,chr(179) f 
,chr(179), 
ichr(179), 
,chr(179), 
,chr(179), 
ichr(179), 
,chr(179), 
,chr(179), 
,chr(I79), 



it 

n 
iV 

6 



1 

Vr 

5 



10 

it 

14 



9 

13 



(★active window*) 

windovC 12,21,80,25); 
end;(*DisplayQAH*) 

(* - -*) 

C* *) 

procedure DiaplayQPSK; 

(*DisplayQPSK shows the QPSK constellation in the upper portion of the 
screen. A window is set at the bottom of the screen for further 
interaction with the user.*) 



begin 

clrscr; 
gotoXY(U f l); 
write(chr(218)); 
for I:= 12. to 64 do 
write(chr(196)); 
write(chr(191)); 
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for I:= 2 to 19 do 
begin 

GdtoXY(ll,I); 
GotoXY(65,I); 
end; 

GotoXY(ll,20); 
write(chr(192)); 
for I := 12 to 64 do 

write(chr(196)); 
write(chr(217)); 
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write(chr(179)); 
writoCchrC 179)); 



1 



window(12,2,63,19); 

writelnC 'This program encodes a QPSK multifrequency signal. ); 
writelnC'The phase ore shown for one frequency.'); 
writelnC.' 
writelnC ' 
writeln( ' 
writelnC 1 
writelnC ! 
writelnC ' 
writelnC ' _ 
writelnC 1 
writelnC 
writelnC ' 
writelnC ' 
writelnC ' 
writelnC 1 
writelnC ' 



* 

3 



frequency. 
' .chrC179),' 
' .cUe(179),' 
\chrC179),' 
',chrC179) ' 
,chrC179) , 
',chrC179), 

.\chrU79), 1 

' chr(179), 
',chrCl79), 
',chrC179),' 
'.chrC^),' 
' ,chrC179),' 
•,chrC179),' 
' ,chrC179),' 



* 

4 



*) 
') 

') 

! 



(^active window*) 

vindowUZ^liSO^); 
end;(*DisplayQPSK*) 

(* *) 

(* - - *) 

procedure SelectQAM; 

(*SelectQAM selects and encodes the symbols for each tone in the band set 
by kl and k2. All symbols enn he randomly or manually selected.*) 

var 

RESPONSE : char; 

ANSWER , I : integer; 

VECTORARRAY. : array[ 0. . 4095] of integer; 

(*- *) 

procedure EncodeData; 

(^EncodeData loads the frequency domain arrays XREAL and XIHAG with the 
QAM symbols in YECTORARRAY and-" their complex conjugate imnr.*.*) 



vnr 



J : integer; 

TEMPR, TEMPI : real; 



begin 

f il lchar (XREALa , s izeof ( XREALa ) , 0 ) ; 
fillchar(XIMAGA,sizeof(XIMAGA),0); 
for J:= kl to k2 do 
begin 

case VECT0RARRAY( J] of 
0 : TEMPI: - 0.0; 
1. : TEMPI:- 90. 0; 
5.. 8 : TEMPI: - 30.0; 



07/29/2004, EAST version: 1.4.1 



31 



5,166,924 



32 



9.. 12 : TEMPI: =-30.0; 

13. . 16: TEMPI: =-90.0; 
end;(*case VECTORARRAY*) 
case VECTORARRAYl J) of 

0 : TEMPR: = 0. 0; 

1,5,9,13 : TEMPR: = 90.0; 

2, 6, 10 ,14: TEMPR: = 30.0; 

3,7,ll,15:TEMPR:=-30.0; 

U t 8, 12, 16: TEMPR: =-90.0; 
ond;(*cosc VECTORARRAY*) 



J ( 



(*Load complex conjugate image") 

XREALa[J]:= TEMPR; 

XBEALa( kx-J] : « TEMPR; . 

XIMAGa(J]:» TEMPI; 

XIMAGAlkx-J]:= -TEMPI; 
end;'(*J:*kl*> 
end; (*EncodeData*) 

(*-— - 

begin 

repeat 
clrscr; 

writeln( 1 Select one of the following for baud 1 , 3AUDC0UNT , ' . 
writeln; 

vriteln('R Randomly select all 1 ,k2 - kl + i,' vectors' 

writeln( 1 1 Individually select vectors 1 ); 

r e ad ln( RESPONSE); 
until RESPONSE in [ 1 r 1 , 1 R' , 1 i' , 1 1' ] ; 
f i 1 lchar ( VECTORARRAY , s izco f ( VECTORARRAY ) , 0 ) ; 
case RESPONSE of 

'r 1 , 'R' : for I:« kl to k2 do 

VECT0RARRAY( I 1 • =r ondom( 1 6 ) + 1 ; 

'i'/l* : begin 

I:« kl; 

while I <" k2 do 
begin 
repeat 
clrscr; 

writcinC Vector Selection Menu'); 
write ln('ff vector desired 1 ); 
writeln('l7 random vector'); 
writeln('lfl no vector 1 ); 
vrite('l9 no more tones Tone ',1,' '); 
readln (ANSWER); 
until ANSWER in [ 1. . 19] ; 
case ANSWER of 

L . 16: VECT0RARRAY( I] : = ANSWER; 

17; VECTORARRAYl I) : 55 random( 16)+1; 
10:;(*nop*) 
19: 1: = k2; 
end;(*case ANSWER*) 
I: =1+1; 

end; (Awhile I*) ? 
end;(*'i' , 'I'*) 
end;(*case RESPONSE 1 *') 
EncodeOata; 
end; (*SelectQAM*) 

(* - - *1 

(• r *j 

procedure SelectQPSK; 
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( ,v SelectQPSK selects end encodes the symbols for each tone in the band 
set by kl* and k2. All symbols can be randomly or manually selected. *) 



var 

RESPONSE : char; 

ANSWER J : integer; 

PHASEARRAY : array[ 0. . 4095] of integer; 

(* : - *) 

procedure EncodeOata; 

(*Encod*Data loads the frequency domain arrays XREAL and XIMAG with the 
QPSK symbols in VECTORARRAY and their complex conjugate image.*) 

var 

J : integer; 

TEHPR, TEMPI : real- 
begin 

f il Ichar ( XREALa , s izeo f ( XREALa ) , 0 ) ; 
fillcharCXIHAGA.sizeofCXIHAGA) ,0); 
for J:= kl to k2 do 
. begin 

case PHASEARRAY[ J] of 
0 : TEMPI: = 0.0; 

1.2 ; TEMPI: = 80.0; 
3,4 i TEMPI: =-80. 0; 

end;(*ea*n PHASEARRAY*) 
case PHASEARRAY[ J] of 

0 :TEHPR:= 0.0; 

1,4 : TEMPR: = 80.0; 

2.3 : TEMPR: =-80. 0; 
end; (*case PHASEARRAY*) 
XREALa( J]:« TEMPR; 
XREALa[ kx-J] : = TEMPR; 
XIMAGa( J]:= TEMPI; 
XIMAGa{ kx-J] : = -TEMPI; 

end;(*J:=kl*) 
end; (*EncodeData*) 

(* - --*) 

begin 
repeat 
clrscr; 

writeln(* Select one of the following for baud f , DAUDCOUNT , 1 . ' ) ; 
writeln; 

writclhC'R Randomly select all 1 ,k2 - kl + l t ' phases'); 

writelnC'l Individually select phases 1 ); 

readln(RESPONSE); 
until RESPONSE in ( 1 r 1 , *R' , ' i 1 , 1 1 ' ) ; 
fillchar( PHASEARRAY, s izcof (PHASEARRAY) ,0); 
case RESPONSE of 

V /R'ifor I;= kl to k2 do 

FHASEARRAY[ I] : =ronriomt + 
'iVl 1 : begin 

I:« kl; 

while I <= k2 do 
begin 

repeat 
clrscr; 

writeln(' Phase Selection Menu'); 

writeln(Js( phase desired* ); 
writeln(*5 random phase 1 ); 



07/29/2004, EAST Version: 1.4.1 



5,166,924 

35 36 

writeln('6 no phase'O; 
write('7 no more tones Tone ' ,2/ '); 
readln (ANSWER); 
until ANSWER in [ 1. . 7] ; 
case ANSWER of 

1. . 4: PHASEARRAY[ I] : = ANSWER; 

5: PHASEARRAY[ I] : = rondom(*)+l; 
6:;(*nop*) 
7: 1: = k2* 
end;(*case ANSWER*) 
I: =1+1; 
end; (*vhile I*) 
end;(* ! i' / I 1 *) 
end;(*case RESPONSE*) 
EncodeData; 
end;(*SelectQPSK*) 

(* - *) 

procedure ScaleData; 

(*ScalcData converts each real valuo in array XREAL down to a byte and 
stores the byte in. packet storage buffer BCST. INDEX sets the location 
in the buffer of each byte*) 



var 

INDEX, J, TEMP : integer; 
DATA :byte; 

begin 

for J :- 0 to kx-1 do 
begin 

TEMP := round( XREALa( J] + 126); 
if TEMP < 0 then 

TEMP := 0; 
DATA : = TEMP; 

INDEX := J.+ (BAUDCOUNT - 1) * kx - 28673; 
BCST[ INDEX) := DATA; 
write(THEFILE , DATA) ; 
end;(*for J*) 
end; (*ScaleData*) 

(* " *) 

(* *) 

procedure Dmainit(vnr DCST : BCSTARRAY; 

BYTECOUNT : integer); 
(^Assembly language procedure used to initialize and unmask the TA\A for 
data transfer. The source must be converted to n BIN file.*) 

external '3MAINIT5. BIN' ; 

(* ;n 

procedure Dmastoo; (*Nasks DMA, stopping data transfer*) 
external 'DMASTCP.BIN'; 

t* *\ 

procedure TINPUT; 

(*This procedure reads in binary file, H0UTFA1. DAT, into an arrav, 
performs. an FFT using CotnpiexFFT and displays the results 
vith GraphData*) 



var 

NUMPTS : integer; 
procodure GetData; * 

(*GetData reads time domain sequence from filo H0UTFA1.DAT into XREAL*) 
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var 

THEFILE : file of byte;. 

DATA :byte; 

I, 

TEMP : integer; 
rdata :real; 

begin 

as sign( THEFILE, 'H0UTFA1. DAT 1 ); 

reset(THEFILE); 

new(XREAL); 

new(XIMAG); 

fillchar(XREALA,sizeof (XREALa),0); 
f il lchar( XIMAGa , s izeo f ( XIHAGa ) , 0) ; 
for I := 0 to kx-1 do 
begin 

read(THEFILE,DATA); (*rend file one byte nt a time*) 
TEMP := DATA; 

temp := temp and 255; (^reduces number of bits used*) 

(*to represent the input data*) 

XREALa[I] := (TEMP - 126); (*londa XREAL army*) 
end; ' (*for I*) 
close(THEFILE); 
end; (*GetData*) 



(* - *) 

(* , - *) 

procedure PutData; 

(''Writes frequency domain arrays XREAL and XINAG to file TFFTOUT. OUT*) 
var 

I : integer; 



OUTFILE ": text; 



begin 

assign(0UTFILE, 'TFFTOUT. DAT' ); 

aewrit e( OUTFILE) 

for I := 0 to kx-1 do 

vritcln( OUTFILE, I, ' * ,XREALa[ I] , 1 1 ,XIMAGa[ I] ); 
ciose( OUTFILE); 
end; (*PutData*) 

( V - 7 - - vV\ 

procedure GraphData; 

(i>GraphData plots the kx/2 points of complex array XREAL nnd XI.MAG. 

Symbols are decoded and represented by the color of each tones 
spectral line, *) 

var 

I, J, (*loop indices") 

KXDIV2, 

COLOR, 

x,xz, - 

Y,YZ : integer; 

begin 

GraphColorMode; 

palette(2); 

Graphbackground( 0) ; 

draw(30, 10,30, 170,1); (*y-axis*) 

draw(30,170,300, 170,1); (*x-axis for full spectrum") 

draw( 30, 155, 300, 155,1); 
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draw(30,80,300,80, 1); (*x-nxis Cor zoom spectrum*) 

(*every other dot equals one bin at 256 Vf ) 
draw(166, 170,166,174,1); (*x-axis scale marks") 
draw(196, 170,196, 17*,1); (*x-axis scale marks*) 
draw(30, 113, 26, 113,1); (*y-exis scale marks*) 
draw(30, 23, 26,23,1); (*y-uxis scale marks*) 

I:= 23; 

(^Grouping of 4 horizontal lines to identify background color 
used for one of the symbols*) 
repeat 

draw(31, I, 300,1,0); 

drav(31, 1+1,300, 1+1,1); ' 

draw(31,I+2,300,I+2,2); 

drawCSlJ+S^OO.I+S^); 

I; =1+12; 
until I > 75; 
KXDIV2 := round(kx / 2); 

(^Decode symbols and assign color*) 
for I : * 0 TO KXDIV2 do 
begin' 

if tXREALA(I]. >= 0) and (XIMAGa( I] > 0) then. 
COLOR := 0; 

if (XREALa(I] < 0) and (XIMAGaU] > 0) then 
COLOR := 1; 

if (XREALa[I] < 0) and (XIHAGa(1] <= 0) then 
COLOR :~ 2; 

if (XREALa[I] >= 0) and (XIHAGa[I] <= 0) then 
COLOR 3; 

i 1 

(*zoom spectrum*) 

if (I >« kl) and .(I <= k2) then 
begin 

XZ ;= rouad((CI * *096.0) / kx) - 1040); 

YZ ;= 80 - round(0.5*(sqrt(sqr(XREALA( I] ) + sqr(XIMAGAi I] )))); 
dmw(XZ , 80 ,X? , YZ , COLOR) ; 
end; 

(*full spectrum*) 

X := round(I/KXDIV2 * 256 +30); 

Y :»170-round(20*ln(sqrt(sqr(XREALA{ I]+i)+ 

sqr( XIHAGa [ I] +l)))/ln( 10)); 
draw(X,170 f X f Y, COLOR); 
end; 

end; (*GrapliData*) 

(* 

begin 

INVERSE := false; (*sets forward FFT*) 
ERROR := 0; 
GetData; 
NUMPTS :* kx; 

write ln( 1 Performing FFT 1 ); 

CoraplexFFT ( NUMPTS , INVERSE , XREAL , XIMAG , ERROR) ; 
writeln(' Error = ' , ERROR, ' hit the enter key' ).; readln; 

PutData; 

GraphData; 
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end; (*Tinput*) 



begin 
repeat 
clrscr; 

assign(TlIEFILE , 1 H0UTFA1. DAT 1 ) ; 

rtvritt(TllEFILE); 

new(XREAL); 

nev(XIHAG); 

INVERSE: =t rue; 

QAM:=false; 

SelectBaud; 

if QAM then • 
"DisplayQAM 
else 

DisplayQPSK; 

(''Packet cont ruction loop*) 

for baudcount : = 1 to numbauds do 
begin 

if QAM then 

SelectQAM 
else 

SelecttjPSK; 
vriteln( 'Performing IFFT* ); 
ComplaxFFT( lex , INVF.RSF, , SURAL . XI MAG t F.RRQR ) ; 
ScaieData; 
end;(*for BAUDCOUNT*) 

dispose(XHEAL); 
dispose(XIMAG); 
close(THEFILE); 

BYTECOUNT := NUMBAUDS*kx-l; ' 

vriteln( ' Press return to transmit 1 ); readln; 

Dmainit ( BCST , BYTECOUNT) ; 

repeat 

writeln( 'Transmit some more? (*yes or no*) 1 ); 

readln( ANSWER); 
until ANSWER in [ V , 'N' , V , 1 Y'] ; 
dmastop; 

window(l,l, 60,25); 
until ANSWER in ( 'n 1 , 'N 1 ] ; 
Tinput; 
end. 
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codcsog segment 

public draainit 
assume cs:codeseg 

procedure DMAINIT ( BLKADDRESS 

DYTECOUNT : 
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DMAINIT 



XMTPOINTER; 
INTEGER); 



44 



this procedure initializes dma channel 3 and sets the 
parameters, to output the array best by passing the address 
of the array start- on the stack. BYTECOUNT is the number 
of bytes to transfer and is pushed on the stack by the 
calling program 



dma equ 
dmapage equ 



0 

80h 



dmainit proc. 
push 
mov 
les 
mov 
out 
out 
mov 
mov 
rol 
push 
and 
out 
pop 
and 
add 
out 
mov 
out 
mov 
out 
mov 
out 
mov 
out 
pop 
ret 

dmainit endp 
codeseg ends 



near 
bp 

J>p,sp 

di t dword ptr( bp+6] 

al,5bh 

dma+ll,al 

dma+12,ni 

ax,es 

d.4 

ax,cl 

ax 

al.Ofh 

dmapage+2,al 
ax 

al.OfOh 

ax,di 

dma+6,al 

al,al\ 

dma+6,el 

ax,(bp+4) 

dma+7 , al 

al ,ah 

dma+7 ,al 

al,3 

dma+l*0,al 

bp 

6 



;use bp to address stack 

;move add of best into es:di 

;dma chan 3 single mode, read, autoinit 

;re«ot first/laat £f 

;caic high order 4 bits of buffer area 



;save ax for dma start addr 
; store in ch 3 dma page rcg 



;get page offset 

; output waveform buffer start addr 



output dma byte count 



; unmask ch 3 to start 



; pop 6 bytes off stack for addr of best 
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APPENDIX £ XMITMES 

program XHITKES; 

(^'Transmits a syncbaud and message from file 'MESSAGE. DAT 1 . The message 
is encoded using QPSK. 'MESSAGE. DAT' is a text file. It should already 
exist before using this program. Output is used to collect data for 
TESTING*) 



const 

FIRSTJLEMENT = -28929; 
type 

TNvector * array[ 0. . 4095] of real; 
TNvectorPtr = ATNvector; 
(*Sync + 61440*) 

BCSTARRAY - array( FIRST_ELEHENT. , 32767] of byte; 

var 

kx, 

kl k2 I w 

NUMBAUDS , HAXNUMB AUDS , 
B AUDCOUNT , BYTECOUNT , 
SYHBOLCOUNT.MAXNUHCHAR, 
HESSAGESIZE : integer; 

MAGNITUDE , 

CHARACTERS_PER_BAUD : real; 
XREAL.XIMAG : TNvectorPtr; 

INVERSE : boolean; 

TEHPDYTE, ERROR : byto; 

BCST : BCSTARRAY; 

BYTEFILE : file of byte; 

TESTFILE : text; 

ANSWER, 

NEXTCltAR :char; 

(*$I FFT87B2. INC*) 
(■*$! COMPFFT. INC*) 

(A......----------- ......... i:] 

procedure SyncDaud; 

(''Process the synchronization baud and stores the 256 point time domain 
sequence at the beginning of the packet storage area.") 

var 

SYNCREAL, SYNCIMAG : TNvectorPtr; 

J, TEMP : integer; 

SYNCDATA : byte; 

SYNCMAG : real; 



begin 

new (SYNCREAL); 
new( SYNCIMAG); 
SYNCMAG: = MAGNITUDE; 

(*load synchronization symbols*) 

f il 1 char( SYNCREALa , s izeo f ( S YNCRE ALa ) , 0 ) ; 

f il lchar( SYNCIMAGa , s izeo f ( SYNCIMAGa ) , 0 ) ; 

SYNCREALa [ 68] : = -SYNCMAG; SYNCIMAGa ( 68] : = -SYNCMAG; 

SYNCREALa( 69} : = -SYNCMAG; SYNCIMAGa ( 69| : = -SYNCMAG; 

SYNCREALa( 70) : = -SYNCMAG; SYNCIHAGa( 70J : = SYNCMAG; 
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SYNCREALa 


711 




= -SYNCMAG; 


SYNCREALa 


72 




= SYNCHAG; 


SYNCREALa 


73 




= SYNCHAG; 


SYNCREALa 


74 




= -SYNCHAG; 


SYNCREALa 


75 




a SYNCHAG; 


SYNCREALa 


76 




- -SYNCMAG; 


SYNCREALa 


77 




= -SYNCHAG; 


SYNCREALa 


78 




■ SYNCHAG; 


SYNCREALa 


79 




.= -SYNCHAG; 


SYNCREALa 


80 




!= SYNCHAG; 


SYNCREALa 


81 




t = SYNCHAG; 


SYNCREALa 


82 




:= -SYNCHAG; 


SYNCREALa 


83 




: «= SYNCHAG; 
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SYNCIMAGa{ 71] : = SYNCMAG 

SYNCIMAGa[72 := -SYNCHAG-, 

SYNCIMAGa( 73] : = SYNCMAG 

SYNCIMAGa[ 74] : = . SYNCMAG 

SYNCIMAGa(75]:= SYtfCMAG 

SYNCIHAGa| 76] : - SYNCMAG 

SYNCIMAGa{ 77] : = -SYNCHAG 

SYNCIMAGa[78]:= -SYNCHAG 

SYNCIHAGa( 79] : = -SYNCMAG 

SYNCIMAGa[ 80] : = SYNCHAG 

SYNCIMAGa[81 := -SYNCMAG 

SYNCIMAGa( 82] : =» • SYNCMAG 

SYNCIHAGa(83 := SYNCMAG; 



(^complex conjugate image*) 
for J :« 68 to 83 do 
begin 

SYNCREALa[ 256-J] : = SYNCREALa{ J] ; 
SYNCIMAGa[ 256-J] : *-SYNCIMAGa[ J] ; 
end;(*for J*) 

CoraplexFFTC 25 6 f INVERSE , SYNCREAL , SYNC I MAG , ERROR) ; 

(*scale data/load time sequence*) 
for J := 0 to 255 do 
begin 

TEMP: "round( SYNCREALa [ J] + 126); 
if TEMP < 0 then 

TEHP:=0; 
SYNCDATA>=TEMP; 

DCST[ J+F1RST_ELEMENTJ : =SYNCDATA; 
end;(*for J*) 
dis pos e( SYNCREAL) ; 
dispose* SYNCIMAG); 
end; (*SyncBaud*) 

c* - 

procedure SelectBaud; 

(*SelectBaud establishes kx, kl, and k2*) 



48 



var 



ANSWER '.integer; 



begin 
kx:=0; 

(*select baud size") 

repeat ' 

if kx < 0 then vriteln('TRY AGAIN'); 

writeln( 'What is the length of the bauds (kx)?'j; 

writein( *i. e. 256, 512, 1024, 2048 % -096".; 

readln( ANSWER); 

case ANSWER of 

256: kx:=256; 

512: kx:=512; 

1024: kx:=1024; 

2048: kx:=2048; 

4096: kx:=4096; 
end; (*case ANSWER*) 
if kx = 0 then kx : = -1; 
until kx > 0; 
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(*set to.ne limits*) 
case kx of 
256: begin 

kl:=68; k2:=83; 
end; 
512: begin 

kl:*135; k2:=166; 
end; 
1024: begin 

kl:=269; k2:=332; 
end; 
2048: begin 

kl:»537; k2:=664; 
. end; 
4096: begin 

kl:=1073; k2:=1328; 
end; 

end; (*case kx ,v ) 
end; (*SelectBaud*) 

(* — 

(* * 

procedure TailorPacket; 

(*TailorPacket sets the raaxiinura number of baud required to encode the 
message*) 

begin 

MESSAGESIZE:* files ize(BYTEFILE); 

vritelnC Message is 1 .MESSAGESIZE, 1 bytes. 1 ); 

(*kx/2 is the number of bit/baud for QPSK. kx/2-2 is the number 

for DQPSK. Each character is 8 bits*) 
CHARACTERS_PER_BAUD:=(kx/8 - 2)/8; 

(*61440/kx is maximum number of bauds possible*) 
MAXNUMCHAR: = trunc( 61440. 0/kx * CHARACTERS PER BAUD); 
if MESSAGESIZE > MAXNUMCHAR then 
begin 

writcln( 'Message is to lnrgn. Tlia last 1 , 
f MESSAGESIZE - MAXNUMCHAR, 

characters will not be transmitted.'); 
MESSAGESIZE: -MAXNUHCHAR; 
end; 

MAXNUMBAUDS; -trunc( MESSAGESIZE / CHARACTERS J'ER.BAUD); 

(^ensure last few characters are included*) 
if frac( MESSAGESIZE / CHARACTERS PER BAUD) > 0. 0 then 
MAXNUMBAUDS: =MAXNUNBAUDS + 1; 

repeat 
uriteln; 

writein( 'Enter number of f ,kx t f bauds to process. ' 

MAXNUMBAUDS , ' is the maximum. 1 ); 
readln(NUMBAUDS); 
until NUMDAUDS in ( 1. . MAXNUMBAUDS] ; 
end; (TailorPacket*) 

(* - 

(* 

procedure Diff Encode; 

(*DiffEncode differential encodes symbols on a tone-to-tone basis. 
BYTEFILE -is read from one byte at -a timo. The byto is isolated into 
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2-bit groups and Atornd in BITS. IUTS is than used to DQPSK onr.ntln th 
frequency domain arrays XREAL and XI HAG. Bytes partially encoded arc 
carried over into the next baud by global variable .TEMPBYTE. '') 



var 



J 

3 ITS 



: integer; 
: byte; 



oegin 

f illchar(XREALA .sizcof (XREALa) ,0) ; 
fillchar(XIHAGA,sizeof (XIMAGa) ,0); 

(*first tone of every baud sec to pi/2*) 
XREALa( kl] : = MAGNITUDE; 
XIMAGa[ kl] : = MAGNITUDE; 

if SYMBOLCOUNT .= 0 then 

read(BYTEFILE , TEHPBYTE) ; 

('''break apart character*) 
for J: = (kl + 1) to k2 do 
begin 

SYMBOLCOUNT: =SYHB0LC0UNT + 1; 

if f rac( SYMBOLCOUNT / 4) = 0.25 then 

BITS:= (TEHPBYTE and $C0) shr 6; 
if f rac( SYMBOLCOUNT / 4) = 0.5 then 

BITS:= (TEHPBYTE and §30) shr A; 
if frac( SYMBOLCOUNT / 4) = 0.75 then 

BITS:= (TEHPBYTE and $0C) shr 2; 
if frac( SYMBOLCOUNT / 4) = 0. 0 then 

begin 

BITS: = TEHPBYTE ond $03; 
if not E0F( BYTEFII.E) then 
read(BYTEFILE .TEHrBYTE) 
else 

TEHPBYTE: =$40; (*fill character*) 
ehd; '• 
if (BITS < 0) and (BITS > 3) then 

writoln( ' Bits not assigned properly'); 
(*differential encode*) 
case BITS of 

0: begin XREALa 



XIMAGa 
1: begin XREALa 
XIMAGa 
2: begin XREALa 
XINAGa 
3: begin XREALa 
XIMAGa 
•nd;(*case BITS*) 
end;(*for J*) 



= XREALa 
- XIMAGa 
=-XIMAGa 
= XREALa 
= XIMAGa 
=-XREALa 
=-XREALa 
=-XIHAGa 



J-l| 
J-l| 
J-l] 
J-l] 
J-l] 
•Ml 
J-l| 
J-1J 



and; 
end; 
end; 
end; 



(^complex conjugate image*) 
for J;= kl to k2 do 
begin 

XREALa [ kx - J]:» XREALa! J]; 
XIHAGa[ kx - J] : =-XIMAGa( J] ; 

writeln(TESTFILE,BAUDCOUNT: 4, J: 5 ,tmnc( XREALa [ J] ): 6, 
trunc(XIMAGA{ J) ):6); 

end; 

end; (*D iff Encode*) 
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(* *) 

(* *) 

procedure ScaleData; 

(*ScaleData converts each real value in XREAL down to a byte ond stores 
the byte in the packet storage buffer, BCST. INDEX establishes the 
location in the buffer of each byte in the packet, *) 

var 

INDEX, J, HHP : integer; 
DATA :byte; 

begin 

for J 0 to kx-1 do 
b«gin 

TEMP := round( XREALa [ J] + 126); 
if TEHP < 0 then 

TEMP := 0; 
DATA := TEMP; 
■ (*256 is added to INDEX to start message bauds 
after the sync baud*) 

INDEX := J+(BAUDCOUNT-l)*kx+FIRST ELEMENT+256; 
BCST[ INDEX] : = DATA; 
end;(*for J*) 
end; (*ScaleData*) 

(* - 

(* *1 

procedure Dmainit(var BCST : BC5TARRAY; 

BYTECOUNT : integer); 

(•'Assembly language procedure used to initialize and unmask the DMA for 

data transfer. The source code must be converted to a BIN file.*) 

external 'DMAINIT5.BIN'; 

(* . Vr\ 

procedure Dmastop; 

(•''Masks DMA, stopping data transfer.*) 
external 'DMASTOP. BIN' ; 

; v - *) 

M> * u) 

begin 

clrscr; 

(^contains hex values to be encoded and transmitted*) 
assign( BYTEFILE , ' MESSAGE. DAT' ) ; 
reset( BYTEFILE) ; 

(^Output file of encoded symbols. Us Ad for syatcm testing*) 
assignCTESTFILE, 1 XMITDAT. DAT' ); 
revrite(TESTTILE); 

INVERSE: =true; 
repeat 

writelnC Enter magnitude of tones, (must be > 0.0) '): 
readln( MAGNITUDE) j 
until MAGNITUDE > 0. 0; 

writelnC 'Loading sync baud.'); 
SyncBaud; 

Select Baud; 
TailorPacket; 
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SYMBOLCOUNT; =0; 
TEMPBYTE:=$00; 

writeln( 'Number of bauds is ' .numbnuds); 

new(XREAL); 

new(XIHAG); 

for baudcount := 1 to nurabauds do 
begin 

DiffEncode; 

writeln( 'Performing IFFT 1 , BAUDCOUNT, r 

NUHBAUDS-BAUDCOUNT, 1 left 1 ); 
ComplexFFT(kx , INVERSE .XREAL, XIMAG , ERROR) ; 
ScaleData; 
end;(*for BAUDC0UNT*) 

dispose(XREAL); 
dispose(XIMAG);. 
close(BYTEFILE); 
close(TESTFILE); 

BYTECOUNT := 256 + NUMBAUDS*kx - 1; 
repeat 

writelnC Press return to transmit' ); rcadln; 

Dmainit( BCST, BYTECOUNT) ; 

repeat 

writeln( 'Transmit 'some more? ( Vf ves or no*) ! ); 

readlnC ANSWER); • 
until ANSWER in [ ! n' ,'N 1 ,'y 1 /Y'J ; 

dmastop; 
until ANSWER in ( V ,'N 1 J; 

(* reset(TESTFILE); 

while not EOF(TESTFILE) do 
begin 

while not EOLN(TESTFILE) do 
begin 

read(TESTFILE ( NEXTCHAR); 
write(NEXTC]IAR); 
end; 

readln(TESTFIL r .); 
writeln; 
end; 

close(TESTFILE);*) 
end. 
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; procedure DMA IN IT ( BLKADDRESS : XHITPOINTER; 

BYTECOUNT : INTEGER); 



;this procedure initializes dma channei 3 and sets the 
;pararaeters to output the array best by passing the address 
5 of the array start, on the stack. DYTECOUNT is the number 
;of bytes to transfer and is pushed on the stack by the 
; calling prograa 

dffla equ 0 
dmapage equ 80h 

dmainit proc. near 

push bp 

mov bp,sp ;use bp to address stack 

les di.dword ptr( bp+6] ; move add of best into es:di 



rol 

push 

and 

out 

pop 

and 

add 

out 

mov 

out 

mov 

out 



mov 
out 
out 
mov 
mov 



al t 5bh 
dma+ll f al 
dma+12, al 
ax,es 
cl,4 
ax,cl 
ax 



;save ax for dma start. addr 



; dma chan 3 single mode, read, autoinit 

; reset first/last ff 

;calc high order 4 bits of buffer area 



mov 
out 
mov 
out 
pop 
ret 



al,0fh 

dmapage+2,al 
ax 

al,0f0h 
ax,di . 
dma+6 ,al 
al,eh 
dma+6 ,al 
ax,(bp+4] 
dma+7 ,al 
al t ah 
dma+7 ,al 
al,3 

dma+l*0,al 

bp 

6 



; get page offset 

'.output uaveform buffer start addr 



; store in ch 3 draa page rcg 



; output dma byte count 



;pop 6 bytes off stack for addr of best 



; unmask ch 3 to start 



dmainit 
codeseg 



endp 
ends 
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APPENDIX C RECMES 

program RECEIVER; 

(''Acquires the signal.. Stores it in a memory buffer. Differential decodes 
between tones. Maximum number of bauds are received. .The number of 
bauds processed is a user input*) 

uses Graph, Crt, tp4d!6; 

(*$!-*) 
(*$R-*) 

const HAXJUFFER = 65500; 
type 

TNvector = array-0. . 4096- of real; (*TYPE for real and imaginary data 

for FFT routing*) 
TNvectorPtr = ATNvector; (*Pointer for FFT data array which 

allows dynamic allocation of memory*) 

var 

INVERSE : boolean; 

XREAL, XIMAG : TNvectorPtr; 
ERROR, TEMPBYTE: byte; 

kl\k2,kx t ANSWER, 
ERR.CODE , BAUDCOUNT, 
SYMBOLCOUNT, 



NUM BAUDS, 
HAXNUHBAUDS : integer; 
MAGNITUDE , 

PHASE : real; 

DATAVECTOR :AintBger; 
DMAPOINTER : pointer; 
TESTFILE, 

OUTFILE :TEXT; 

FFT87B2. INC*) 
COHPFTT. INC*) 

(* *) 

(* - ' *) 



procedure PacketSetUp; 

begin 
repeat 
clrscr; 

if kx < 0 then writeln( f TRY AGAIN 1 ); 
writeln( 'Enter baud size '); 
readlnC ANSWER); 
case ANSWER of 
256: kx: = 256; 

512: 1 kx: = 512; : 
1024: kx: =1024; 
2048: kx: =2048; * 
4096: kx: =4096; 
end; (*case*) 

if kx = 0 then kx := -1; 
until kx > 0; 

HAXNUHBAUDS := trunc((MAXJUJFFER/2)/kx); 
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repeat 
writelri; 

writeln( 'Enter number of '.kx, 1 bauds to process. 

MAXNUMBAUDS , * is the maximum.'); 
readln(NUMBAUDS); 
until NUMBAUDS in ( 1. . HAXNUMBAUDSj ; 

ki ;-,round(kx * 67. 0 / 256.0 + 1); 
k2 i— round(kx* 83.0 / 256.0); 
end; (*PaeketSetUp*) 

(*-— 

(* - 

procedure AcquireData; 

(*AcquireData initializes Metrobyte DASH-16F data acquisition board, 
usin* TT00LS procedure D16_int and D16_ainm. Data transfer is 
controlled by the DMA controller and initialized by D16_ainm and 
disabled by D16_draa_int disable. TT00LS procedures arc external 
procedures included by uses' tp4dl6.*) 



var 

RATE: real; 

I,CtfT_NUM, MODE , CYCLE, TRIGGER, 
BASE.AOR, INT.LEVEL, DMA_LEVEL, 
BOARD JWM, CHANLO, 

OPJTYPE, STATUS , NEXT.CNT, ERR_C0DE_S : integer; 
begin 

BOARD.NUH := 0; INT.LEVEL := 7; DMAJLEVEL := 1; 
BASE_ADR := $300; 

D 16_init ( BOARD J<UM , B ASE_ADR , INT.LEVEL , DMA.LEVEL , ERR.CODE ) ; 



CHANLO := 0; 

CYCLE: "0; (*0-one sweep of the DMA 1-autoinitializc*)' 
TRIGGER: *0; (*0 - external'* 1 - internal*) 
CNT_NUM:=trunc(HAX - BUFFER / 2); of samples*) 
RATE :■ 10000- 0; (*used for internal trigger*) 
MODE := 2; (*DMA mode* ) 
writeln( 'Ready to acquire ); 

D16_ainm( BOARDJWH , CHANLO .MODE .CYCLE .TRIGGER ,CNT_NUM , 
RATE , DATAVECTORa , ERR_COUE ) ; . 

STATUS := 11; 

(*status indicates the progress of acquisition. When nil 
samples have been acquired stotus c 0*) ( 
repeat ! 

D16_dma JLnt_s t atus( BOARD JWN ,or_TYPE , STATES , NEXTJINT, K RK _C0DE_S) ; 
until STATUS « O; 
writeln( 'Data received 1 ); 

if ERRJ30DE <> 0 then 

D16_print_error( ERR.CODE) ; 
D16_dma_int_disable( B0ARD_NUM , ERR_C0DE) ; 

end; (* Acquire*) 



procedure ConvertData; 

(*ConvertData separates channel and acquired data. CHAN.DATA is not used. 
Acquired data is stored in XREAL. *) " 
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var 

AD.DATA: array( 0. . 4096] of integer; 
I , CHAN.DATA , ERR_C0DE , 
SEGMENTPART , OFFSETPART: integer; 
NEWDATAVECTOR : Ainteger; 

TEMPPOINTER : pointer; 

begin 

writeln( 'Processing baud ' .BAUUCOUNT); 
fillchar(XREALA,si2eof (XREALa),0); 
£illchar(XIMAGA,si2eaf (XIMAGa),0); 
SEGMENTPART: =s eg( DATAVECTORa ) ; 

OFFSETPART: *=ofs(DATAVECTOR a) + 2 * kx * (BAUDCOUNT - 1); 
TEMPPOINTER: =pt r ( SEGMENTPART f OFFSETPART) ; 
NEWDATAVECTOR : « TEMPPOINTER; 

dl6_convert_dat a( 2047 , kx , NEWDATAVECTORa , AD J)ATA[ 0] , 
CHAN_DATA , 0 , ERR_CODE ) ; . 

for I: » 0 to (kx - 1) do * 
begin 

XREALa[ I] : e ADJ)ATA[ I] /5; 
end; 

end; (*ConvertData*) 

(* »«- ft) 

procedure DiffDecode; 

<*DiffDecode dif f erentially decodes complex frequency domain arrays XREAL 
and XIMAG. Four decoded symbols ere recombined into a byte -and 
transferred to file BYTES OUT. DAT. *) 

var 

I : integer; 

TEMPRE AL , TEMP I MAG : real; 
BITS :byte; 
TEMPCHAR : char; 



begin 

for I := kl to (k2 - 1) do 

begin . ... 

(^Complex multiply two adjaeoiit tonon , r and t:ho complex cmijngnt.n of Ml. 
Tills will givo the phase difference betwncn chn two tones, Th* .mswo.r 
is in rectangular notation*) 

TEMPREAL:»XREALa( 1] * XREALa[ + XIMAG a( I] * XIMAGaI 

TEMPIMAG:=XREALa[I] * XIMAGa( 1+1] - XREALa[ I+l| * XIMAGa( I| ; 

(^Complex multiply (TEHPREAL + j TEMP IMAGE) and (1+j). This rotates the 
differential vector pi/A radians, XREAL and XIMAG are used to score 
the results. This eliminate the original data*) 

XREALa[I]:=(TEMPREAL - TEMPIMAG)/80; 

XIHAGa( I] : =( TEHPREAL + TEMPIMAG)/80; 

(^decode*) 

if (XREALa(I) >= 0) and (XIMAGaI I] > 0) then 
BITS: =$00; 

if .(XREALa(I] < 0) and (XIMAGa(I1 > 0) then 
BITS: =$01; 

if (XREALa(I) < 0) and (XIMAGa[I] <= 0) then 
BITS: =$03; 

if (XREALa[I) >= 0) and (XIMAGa( IJ <= 0) then 
BITS: =$02; 

SYMB0LC0UNT : = SYMBOLCOUNT + 1; 
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(*£I11 TEMPBYTE with four symbals*) 
if f rac( SYMBOLCOUNT / 4) » 0.25 then 

TEMPBYTE := (BITS shl 6); 
if frac( SYMBOLCOUNT / 4) = 0.5 then 

TEMPBYTE := (BITS shl 4) or TEMPBYTE; 
if f rac( SYMBOLCOUNT / 4) = 0.75 then 
;.. TEMPBYTE (BITS shl 2) or TEMPBYTE; 
if ( frac( SYMBOLCOUNT / 4) * 0: 0) then 

begin 

TEMPBYTE := BITS or TEMPBYTE; 
TEMTCIIAR : => chr( TEMPBYTE); 
write( OUTFILE , TtHPCHAR) ; 
TEMPBYTE: =0; 
end; (*if frac*) 
end; (*for I*) 

XREALA[k2J ;=l. 0; 
XIMAGa(1c2):=1.0; 

end; (*Diff Decode*) 



procedure Showmessage; 

(*Showraessage reads in decoded message*) 

Var NEXTCIlAR: char; 



begin 

writeln; 

writeln('The message transmitted is.. 1 ); 

assign( OUTFILE, 'MESSAGE. DAT 1 ); 

rn«*t(OUTFII,E); 

while not E0F( OUTFILE) do 

begin * ■ 

while not EOLN( OUTFILE) do 
begin 

read( OUTFILE , NEXTCIIAR) ; 

writ e( NEXTCHAR) ; 
end;(*while not EOLN*) 
readln( OUTFILE); 
writeln; 

•nd;(*while not EOF*) 
ciose(OUTFILE); 
end; (*Showmessage*) 

(* - - *) 

begin (*main body*) 

Ge tOMABuf f er( HAX_BUFFER , DMAPOINTER , ERR_C0DE ) ; 

DATAVECTOR := DMAPOINTER; (*This statement assigns a generic pointer 

to a variable of a specific pointer typo, 
i. e. -Ainteger, so that the pointer cnn be 
passed to the d!6_ainm routine.*) 

new(XREAL); 
new(XIMAG); 



INVERSE := false; (*snts forward FFT*) 

ERROR 0; 
kx: =0; 
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PacketSetUp; 

(♦received message output file ,v ) 
ass ign( OUTFILE , ' MESSAGE. DAT* ) ; 
rewrite(OUTFILE); 

(*RECSTAT file shows the computed real and Imaginary values of the 
received signal and its decoded representation''') 

assign(TESTFILE,' RECSTAT.DAT'); 

rewrita(TESTTILE); 
(* writeln(TESTFILE,'Baudlength is' ,kx:S);*) 

SYMBOLCaUNT:=0; 
TEHPBYTE: =0; 

AcquireOata; (*AcquiroData samples input analog signal'') 

for BAUDCOUNT := 1 to NUMBAUDS do 
begin 

ConvextOata; 

ComplexFFT (kx, INVERSE, XREAL.XIHAG. ERROR); 
(* writelnCTESTFILE.'Baud': 11. 'Tone': 5. 'Real' : 10, 

1 Imaginary ' : 10 , ' Mag ' : -10 , ' Phase ' : 8) ; *) 
for J:=kl to k2 do 
begin 

MAGNITUDE: =20* ln(sqrt(sqr(XREAI,A( J] )+ 
( sqr(XIHAGA[ JJ ))) / ln(10); 

if (XREALa{J1 >= 0) then 

PHASE: =arctan(XIHAG A [ J) / XREALa(JI); 

if (XREALa(J) < 0) and (XIMAGa[J] > 0) then 
PHASE: =Pi + arctan(XIHAGA( J] / XR£ALa( J] ); 

if (XREALa(J) < 0) and (XIHAGa(J) <= 0) then 
PHASE: =arctan(XIMAGA[JJ / XREALa(JJ) - Pi; 

writeln(TESTFILE,'R' , BAUDCOUNT: 4,J: 5 , 

XREALa( J] : 10: 4,XIMAGa( J] : 10: 4 , 
MAGNITUDE: 10: 4,ri!ASE: 8: 4)", 

end; 

Diff Decode; 

for J:«kl to k2 do 
begin 

MAGNITUDE: -sqrt(8qr(XREALA( J] )-faqr(XIHAGA{ JJ )); 

writeln(TESTFILE /D 1 , BAUDCOUNT: 4 , J:-5 ; 

XREALa[ J] /MAGNITUDE: 10: 4, 
XIMAGa[ J] /MAGNITUDE: 10: 4); 

end; 

end; 

writeln( OUTFILE , chx(TEMPBYTE) ) ; 

close(OUTFILE); 
close(TESTFILE); 
dispoi.(XREAL); 
dispose(XIHAG); 

FreeDMABuff er( MAX_BUFFER , DMAPOINTER , ERR_C0DE) ; 
Showmesaage; 

(* writeln( Error = * , ERROR,' hit the enter key'); read In;*) 
end. 



07/29/2004, EAST Version: 1.4.1 



5,166, 

69 

I claim: 

1. A data-driven, two-stage echo cancellation circuit 
for use in a full duplex, communications system which 
encodes and decodes , transmission signal data in the 
time domain and the frequency domain, said echo can- 5 
cellation circuit comprising: 

a first stage for generating a first estimate of an echo 
signal utilizing transmission signal data in the time 
domain, said first estimate of an echo signal sub- 
tracted from a received analog signal; and 10 

a second stage, for generating a second estimate of 
residua] echo signal utilizing transmission signal 
data in the frequency domain, said second estimate 
subtracted from said received signal after said re- 
ceived signal has been converted to a digital signal 15 
in the frequency domain. 

2. An echo cancellation circuit as in claim 1 wherein 
said first stage comprises: 

a linear filter having an input for receiving said trans- 
mission signal data, said linear filter generating said 20 
first estimate of an echo signal in digital form; 

a digital-to-analog converter coupled to an output of 
said linear filter for converting said digital first 
estimate to an analog form; and 

a summer circuit coupled to an output of said digital- 
to-analog converter and to said communications 
channel, said summer circuit for subtracting said 
first estimate from said received analog signal. 

3. An echo cancellation circuit as in claim 2 wherein 
said second stage comprises: 

a multiplier circuit having first and second inputs, 
said first input for receiving said transmission signal 
in the frequency domain, said second input for 
receiving complex weighting coefficients, said 35 
transmission signal multiplied with said weighting 



,924 

70 

coefficients for generating said second estimate of 
said residual echo signal; and 
a summer circuit coupled to an output of said multi- 
plier circuit and to said communications channel, 
said summer circuit for subtracting said second 
estimate of residual echo signal from said digital 
received signal in the frequency domain. 

4. An echo cancellation circuit as in claim 2 wherein 
said linear filter comprises a finite impulse response 
linear filter implementing a weighting function, h(i), 
defined by the following equation: 

A fl +J(0=A fl (0+* J r- 1 (K«Wfl-0) 

i=0,l...N-l, n=0,l ...k*-l 
where i-the number of filter coefficients, N=one of a 
number of filter coefficients, n=the number of sample 
points in a baud of said transmission signal, k x =one of 
a number of digital signal samples, r(n)=the receive 
time domain digital signal sequence and x(n)=the trans- 
mit time domain digital signal sequence. 

5. An echo cancellation circuit as in claim 3 wherein 
said complex weighting coefficients, C(k), are defined 
by the following equation: 

a*) = i Rmxm/ s \xm\ . 

k s k\ . . . *2. where k » the number of 

sample points in a baud of said transmission signal, 
Xi(k)= values in the transmit frequency domain array 
and Ri(k)= values in the receive frequency domain 
array. 

***** 
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